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Preface 

A textbook presents far more material than any professor can cover in class. 
These lecture notes present only somewhat more than I covered during the 
half-semester course Spectral Theory of Partial Differential Equations (Math 
595 STP) at the University of Illinois, Urbana-Champaign, in Fall 2011. 

I make no claims of originality for the material presented other than some 
originality of emphasis: I emphasize computable examples before developing 
the general theory. This approach leads to occasional redundancy, and some- 
times we use ideas before they are properly defined, but I think students gain 
a better understanding of the purpose of a theory after they are first well 
grounded in specific examples. 

Please email me with corrections, and suggested improvements. 



Richard S. Laugesen Email: Laugesen@illinois.edu 

Department of Mathematics 

University of Illinois at Urbana-Champaign, U.S.A. 



Prerequisites and notation 

We assume familiarity with elementary Hilbert space theory: inner product, 
norm, Cauchy-Schwarz, orthogonal complement, Riesz Representation The- 
orem, orthonormal basis (ONB), bounded operators, and compact operators. 
Our treatment of discrete spectra builds on the spectral theorem for compact, 
selfadjoint operators. 

All functions are assumed to be measurable. We use the function spaces 

L = integrable functions, 
L^ = square integrable functions, 
L°° = bounded functions, 

but we have no need of general L'^ spaces. 

Sometimes we employ the L^-theory of the Fourier transform. 



f(y 



flxje-^^^^-^'dx. 



Only the basic facts are needed, such as that the Fourier transform preserves 
the L^ norm and maps derivatives in the spatial domain to multipliers in the 
frequency domain. 

We use the language of Sobolev spaces throughout. Readers unfamiliar 
with this language can proceed unharmed: we mainly need that 

H = W ' = {L -functions with 1 derivative in L }, 
H = W ' = {L -functions with 2 derivatives in L }, 

and 

HJ = Wq'^ = {H^ -functions that equal zero on the boundary}. 

(These characterizations are not mathematically precise, but they are good 
enough for our purposes.) Later we will recall the standard inner products 
that make these spaces into Hilbert spaces. 

For more on Sobolev space theory, and related concepts of weak solutions 
and elliptic regularity, see |Evansj . 



Introduction 

Spectral methods permeate the theory of partial differential equations. One 
solves linear PDEs by separation of variables, getting eigenvalues when the 
spectrum is discrete and continuous spectrum when it is not. Linearized 
stability of a steady state or traveling wave of a nonlinear PDE depends on 
the sign of the first eigenvalue, or on the location of the continuous spectrum 
in the complex plane. 

This minicourse aims at highlights of spectral theory for selfadjoint par- 
tial differential operators, with a heavy emphasis on problems with discrete 
spectrum. 

Style of the course. Research work differs from standard course work. 
Research often starts with questions motivated by analogy, or by trying to 
generalize special cases. Normally we find answers in a nonlinear fashion, 
slowly developing a coherent theory by linking up and extending our scraps 
of known information. We cannot predict what we will need to know in 
order to succeed, and we certainly do not have enough time to study all 
relevant background material. To succeed in research, we must develop a 
rough mental map of the surrounding mathematical landscape, so that we 
know the key concepts and canonical examples (without necessarily knowing 
the proofs). Then when we need to learn more about a topic, we know where 
to begin. 

This course aims to develop your mental map of spectral theory in partial 
differential equations. We will emphasize computable examples, and will be 
neither complete in our coverage nor completely rigorous in our approach. 
Yet you will finish the course having a much better appreciation of the main 
issues and techniques in the subject. 

Closing thoughts. If the course were longer, then we could treat topics 
such as nodal patterns, geometric bounds for the first eigenvalue and the 
spectral gap, majorization techniques (passing from eigenvalue sums to spec- 
tral zeta functions and heat traces), and inverse spectral problems. And we 
could investigate more deeply the spectral and scattering theory of opera- 
tors with continuous spectrum, giving applications to stability of traveling 
waves and similarity solutions. These fascinating topics must await another 
course. . . 
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Part I 



Discrete Spectrum 



Chapter 1 
ODE preview 



Goal 

To review the role of eigenvalues and eigenvectors in solving 1st and 2nd order 
systems of linear ODEs; to interpret eigenvalues as decay rates, frequencies, 
and stability indices; and to observe formal analogies with PDEs. 



Notational convention 

Eigenvalues are written with multiplicity, and are listed in increasing order 
(when real-valued) : 

Ai < A2 < A3 < ■ ■ • 

Spectrum of a real symmetric matrix 

If A is a real symmetric d x d matrix {e.g. A = [g^] when d = 2) or 
Hermitian matrix then its spectrum is the collection of eigenvalues: 

spec(A) ={Ai,...,Ad}c M 

(see the figure). Recall that 

Avj = AjVj 

where the eigenvectors {vi, . . . ,Vd} can be chosen to form an ONB for M'^. 



Ai A2 
spec(A) 



C 



^Xr 



^Xr 



Observe A : M'^ ^ M'^ is diagonal with respect to the eigenbasis: 







Cl 



Cd 



AiC 



l"-! 



AdCd 



What does the spectrum tell us about linear ODEs? 

Example 1.1 (1st order). The equation 

— = -Av 
dt 

v(0) = Y_ CjVj 
has solution 

Notice Aj =decay rate of the solution in direction Vj if Aj > 0, or growth 
rate (if Aj < 0). 

Long-time behavior: the solution is dominated by the first mode, with 

v(t) ~ e^^^^CiVi for large t, 

assuming Ai < A2 (so that the second mode decays faster than the first). The 
rate of collapse onto the first mode is governed by the spectral gap A2 — Ai 



9 



since 

d 

~e-^'*(ciVi+0(e-*^2-Ai)t^_ 

Example 1.2 (2nd order). Assume Ai > 0, so that all the eigenvalues are 
positive. Then 

—J =-Av 
vfO) = > CjVj 






has solution 



v(t) = cos(v^t)v(0) + —= sm{y/At)V{Q) 

VA 



= y cos(A/AJt)CjVj + 2~ ^^sin(A/AJt)c|vj. 

V ^ 

Notice ^yXj =frequency of the solution in direction Vj. 

Example 1.3 (1st order imaginary). The equation 

.dv , 
I—- = Av 
dt 



V|U) = > CjVj 



has complex-valued solution 

v(t) = e-^'^V(0) = _^e-^''Vj- 
This time Aj =frequency of the solution in direction Vj 

10 



What does the spectrum tell us about nonlinear ODEs? 
In/stability! 

Example 1.4 (1st order nonlinear). Suppose 



dv 
dt 



Ffv) 



where the vector field F satisfies F(0) = 0, with first order Taylor expansion 

F(v) =Bv + 0{\vf] 

for some matrix B having d linearly independent eigenvectors Vi , . . . , v^ and 
corresponding eigenvalues Ai, . . . , A^ G C (The eigenvalues come in complex 
conjugate pairs, since B is real.) 

Clearly v(t) = is an equilibrium solution. Is it stable? To investigate, 
we linearize the ODE around the equilibrium to get ^ — Bv, which has 
solution 

v(t) = e^\iO) = ^e^''cjVj. 

Notice v(t) ^ as t ^ oo if Re(Aj) < for all j, whereas |v(t)| ^ oo if 
Re(Aj) > for some j (provided the corresponding coefficient Cj is nonzero, 
and so on). Hence the equilibrium solution v(t) = is: 

• linearly asymptotically stable if spec(B) C LHP, 



X 



^< X^ 



spec(B) 



C 



• linearly unstable if spec(B) fl RHP ^ 0. 

The Linearization Theorem guarantees that the nonlinear ODE indeed be- 
haves like the linearized ODE near the equilibrium solution, in the stable 
and unstable cases. 
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c 



^<r 



^<r 



spec(B) 



The nonlinear ODE's behavior requires further investigation in the neu- 
trally stable case where the spectrum lies in the closed left half plane and 
intersects the imaginary axis (Re(Aj) < for all j and Re(Aj) = for some 



^< >^ 



C 



spec(B) 



For example, if B = [^ q^] (which has eigenvalues ±1), then the phase 
portrait for ^ = Bv consists of circles centered at the origin, but the phase 



portrait for the nonlinear system ^ = Y[v] might spiral in towards the origin 
(stability) or out towards infinity (instability), or could display even more 
complicated behavior. 



Looking ahead to PDEs 

Now suppose A is an elliptic operator on a domain D. C M'^. For simplic- 
ity, take A = —A. Assume boundary conditions that make the operator 
self-adjoint (we will say more about boundary conditions later). Then the 
eigenvalues Aj and eigenfunctions Vj (x) of the Laplacian satisfy 



-Avj = AjVj 



in a 



and the spectrum increases to infinity: 

Ai < A2 < A, < ■ ■ ■ 



00. 
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The eigenfunctions form an ONB for L (II). 

Substituting A = —A into the ODE Examples 1J_ 1^ transforms them 
into famous partial differential equations for the function v(x, t). We solve 
these PDEs formally by separation of variables: 



Example 1.1 — diffusion equation Vt = Av. Separation of variables 



gives the solution 

v = e%(-,0) = }^e-^iVj 
where the initial value is v(-,0) = Y. CjVj. Here Aj =decay rate. 



Example 1.2 — wave equation Vtt = Av. Separation of variables gives 



V = cos( V^t)v(-, 0) + -^^ sin( V^t)vt(-, 0) 

= 2_ cos(A/AJ^t)CjVj + 2_ —j= sinfA/AJ'tjcjVj. 

V J 

So a/Aj^ =frequency and Vj =mode of vibration. 



Example |1.3| — Schrodinger equation ivt = — Av. Separation of vari- 
ables gives 

Here Aj ^frequency or energy level, and Vj =quantum state. 



We aim in what follows to analyze not just the Laplacian, but a whole 
family of related operators including: 



= -A 


Laplacian, 


= -A + V(x) 


Schrodinger operator 


= (iV + A)^ 


magnetic Laplacian, 


= (-A)^ = AA 


biLaplace operator. 



The spectral theory of these operators helps explain the stability of different 
kinds of "equilibria" for evolution equations: steady states, standing waves, 
traveling waves, and similarity solutions. 
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Chapter 2 

Laplacian — computable 
spectra 

Goal 

To develop a library of explicitly computable spectra, which we use later to 
motivate and understand the general theory. The examples are classical and 
so proofs are left to the reader, or else omitted, except that Weyl's asymptotic 
law is proved in detail for rectangles. 

References [Strauss] Chapters 4, 10; [Farlow] Lesson 30. 

Notation 

Let O be a bounded domain in M'^, d > 1 . Fix L > 0. 
Abbreviate "boundary condition" as "BC" : 

• Dirichlet BC means u = on 90, 

• Robin BC means |^ + cru = on 90 (where a G M is the Robin 
constant), 

. Neumann BC means f^ = on 90. 

Spectra of the Laplacian 

/ 9 \2 ' ' 

A = V-V= (:— ) +■■■ + 



\9xi/ v9xd 
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Eigenfunctions satisfy — Au = Au , and we order the eigenvalues in increas- 
ing order as 

Ai < A2 < A3 < ■ • 



00. 



To get an ONB one should normalize the eigenfunctions in L^, but for sim- 
plicity, we will not normalize the following examples. 



One dimension —u" = Au 



1. Circle T = R/lnZ, periodic BC: u(-7r) = u(7r),u'(-7t) = u'in). 

Eigenfunctions e''^" for j G Z, or equivalently 1, cos(jx), sin(jx) for j > 1. 
Eigenvalues Aj = f for j G Z, or A = 0^ V,V, 2\ 2^ . . . 



2. Interval (0,L) 

(a) Dirichlet BC: u(0) = u(L) = 0. 

Eigenfunctions iij(x) = sin(J7rx/L) for j > 1. 
Eigenvalues Aj = (J7t/L)^ for j > 1, e.g. L = n 



1^ 2^ 32 



Dirichlet spectrum 



^<r 



^<r 



^<r 



2/_2 



AL^TT 




(b) Robin BC: -u'(0) + 0u(O) = u'(L) + ctu(L) = 0. 
Eigenfunctions Uj(x) — y^cos(y^x) + 0sin(y^x). 

Eigenvalues Pj — jth positive root of tan(y^L) = -^^ for j > 1. 

(c) Neumann BC: u'(0) = u'(L) = 0. 

Eigenfunctions Uj(x) = cos(J7rx/L) for j > (note uq = 1). 
Eigenvalues [i^ = [jn/L]^ for j > 0, e.g. L = n =^ A = 0^ ^^2^3^ 
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Robin spectrum 



pLV 




Neumann spectrum 



>^i«^ 



lA?l7i^ 




Spectral features in 1 dim 

i. Scaling: eigenvalue must balance d^/dx^, and so A ~ (length scale)^^. 

Precisely, Aj((0,tL))=Aj((0,L))/t2. 
ii. Asymptotic: eigenvalues grow at a regular rate, Aj - (const. )j^ 
iii. Robin spectrum lies between Neumann and Dirichlet: 

Neumann < Robin > Dirichlet 

as one sees formally by letting a approach or oo in the Robin BC |^ + cru = 
0. 



Two dimensions 



-Au = Au 



1. Rectangle Q = (0, L) x (0, M.) (product of intervals). 

Separate variables using rectangular coordinates Xi,X2. See the figures at 
the end of the chapter! 

(Note that every rectangle can be reduced to a rectangle with sides par- 
allel to the coordinate axes because the Laplacian, and hence its spectrum, 
is rotationally and translationally invariant.) 

(a) Dirichlet BC: u = 

Eigenfunctions Ujic(x) — sin(J7rxi/L) sin(k7tX2/M) for j,k > 1. 
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Eigenvalues Aji^ = []n/L)^ + [kn/M]^ for j,k > 1, 
e.g.L = M = n =^ A = 2,5,5, 8, 10, 10, . . . 

(b) Neumann BC: |^ = 

Eigenfunctions Ujk(x) = cos(J7txi/L) cos(k7tX2/M) for j,k > 0. 
Eigenvalues \x^ = [jn/L]^ + [kn/M]^ for j,k > 0, 
e.g.L = M = n =^ A = 0, 1, 1,2,4,4, .. . 

2. Disk a = {x G M^ : |x| < R}. 

Separate variables using polar coordinates r, 9. 

(a) Dirichlet BC: u = 
Eigenfunctions 

Jo(rjo,m/l^) for m> 1, 

]niT]n,mf/^] cos[nQ] and jTi(TJn,mVl^) sin(n9) for n > 1,tu > 1. 

Notice the modes with n = are purely radial, whereas when n > 1 the 
modes have angular dependence. 

Eigenvalues A = (jn,m/R)^ for n > 0, ra > 1 , where 

Jn = Bessel function of order n, and 
jn,m = TR-th positive root of Jn(T) = 0. 

The eigenvalue An^in has multiplicity 2 when n > 1 , associated to both cosine 
and sine modes. 

From the graphs of the Bessel functions Jo, Ji, Ji we can read off the first 
4 roots: 

jo,i ~ 2.40, ji,i ~ 3.83 J2,i ~ 5.13 ji,2 ~ 5.52. 

These roots generate the first 6 eigenvalues (remembering the eigenvalues are 
double when n > 1 ) . 

(b) Neumann BC: |^ = 

Use roots of j;(r) = 0. See [Bandlel Chapter III]. 

3. Equilateral triangle of sidelength L. 

Separation of variables fails, but one may reflect repeatedly to a hexagonal 
lattice whose eigenfunctions are trigonometric. 

Dirichlet eigenvalues Aji^ = ^^ [f + jk + k^) for j, k > 1 . 

Neumann eigenvalues |J^k = ^^i]^ + jk + k^) for j, k > 0. 
See [Mathews fc Walker j [Mc( Martin] . 
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Spectral features in 2 dim 

i. Scaling: eigenvalue must balance A, and so A ~ (length scale)^^. 
Precisely, Aj(ta) = Aj(a)/t^. 

ii. Dirichlet and Neumann spectra behave quite differently when the 
domain degenerates. Consider the rectangle, for example. Fix one side length 
L, and let the other side length M tend to 0. Then the first positive Dirichlet 
eigenvalue blows up: taking j = k = 1 gives eigenvalue [n/L]^ + (tt/M)^ -^ 
oo. The first positive Neumann eigenvalue is constant (independent of M): 
taking j = 1 , k = 0, gives eigenvalue (tt/L)^. 

iii. Asymptotic: eigenvalues of the rectangle grow at a regular rate. 

Proposition 2.1. (Weyl's law for rectangles) The rectangle (0, L) x (0, M.) 
has 

-\.^yi^^ as 3 ^ OO, 

Area 

where Area — LM is the area of the rectangle and Ai , A2, A3, . . . and [Xi , 1^2, IJ.3, . , 
are the Dirichlet and Neumann eigenvalues respectively, in increasing order. 

Proof. We give the proof for Dirichlet eigenvalues. (The Neumann case is 
similar.) Define for a > the eigenvalue counting function 



N ( (x) = ^{eigenvalues < ex} 



r 



#{(j,k)GNxN:(j,k)GE} 



<!} 



where E is the ellipse (x/a)^ + (y/b)^ < 1 and a = -/ah/yt, b = y/RM/n. 




E-(l,l) 



X X X X X 



X X X X X 



X X X X X 




We associate each lattice point (j,k) G E with the square 

S(j,k) = [j-l,j]x [k-1,k] 

whose upper right corner hes at (j,k). These squares all lie within E, and so 
by comparing areas we find 

N((x) < (area of E in first quadrant) = -Trab = a. 

4 4n 

On the other hand, a little thought shows that the union of the squares covers 
a copy of E shifted down and left by one unit: 

U(j,k)eES(j,k) D (E— (1, 1)) n (first quadrant). 

Comparing areas shows that 

N ( a) > T^tab — a — b 
4 

LM L+M ^ 
-oc vo. 



47T n 

Area Perimeter ^- 

Combining our upper and lower estimates shows that 

^ , , , Area 
N((x)~^^a 

as a ^ oo. To complete the proof we simply invert this last asymptotic, 
with the help of the lemma below. D 



Lemma 2.2. (Inversion of asymptotics) Fix c > 0. Then: 

a, 
c 



Nice) 



Aj~cj. 
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Proof. Formally substituting a = Aj and Nfoc) = j takes us from the first 
asymptotic to the second. The difficulty with making this substitution rig- 
orous is that if Aj is a multiple eigenvalue, then N(Aj) can exceed j. 

To circumvent the problem, we argue as follows. Given e > we know 
from N[a) ~ cx/c that 

(1 -£)-<N(a) <(!+£)- 
c c 

for all large oc. Substituting a = Aj into the right hand inequality implies 
that 

c 
for all large j. Substituting a = Aj — 5 into the left hand inequality implies 
that 

c 
for each large j and < 5 < 1 , and hence (by letting 8 ^> 0) that 

(l-e)^<j. 
c 



We conclude that 



<4< 



1 + £ CJ 1 — £ 

for all large j, so that 

lim -^ — 1 

j->oo C) 

as desired. D 



Later, in Chapter 11, we will prove Weyl's Asymptotic Law that 

Aj ~ Anj/Aiea. 

for all bounded domains in 2 dimensions, regardless of shape or boundary 
conditions. 

Question to ask yourself What does a "typical" eigenfunction look like, 
in each of the examples above? See the following figures. 
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Dirichlet mode (j=2,k=3) 



Dirichlet square (L=M=7r) 




Neumann mode (j=2,k=3) 



Neumann square (L=M=7r) 

w ^ ^ ^ 




Dirichlet sine mode (n=2,m=l) 



pirichlet disk (R=^fn ,A=7i^] 

X n n X n 
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Chapter 3 

Schrodinger — computable 
spectra 

Goal 

To study the classic examples of the harmonic oscillator (1 dim) and hydrogen 
atom (3 dim). 



References [Strauss] Sections 9.4, 9.5, 10.7; GustafsonSigal Section 7.5, 
7.7 



Harmonic oscillator in 1 dimension 



— u" + X u = Eu 



Boundary condition: u(x) ^ as x -^ ±oo. (Later we give a deeper per- 
spective, in terms of a weighted L^-space.) 

Eigenfunctions Ui<(x) = H\^[x)e^^ ^'^ for k > 0, where Hi^ = k-th Hermite 
polynomial. 

Eigenvalues E^ = 2k + 1 for k > 0, or E = 1 , 3, 5, 7, . . . 
Examples. Ho(x) = l,Hi(x) =2x,H2(x) = 4x^-2, Hk(x) = (-Ij'^e"' (|^) V"' 

Ground state: Uo(x) = e^" ^^ = Gaussian. (Check: — Uq + x^Uq = Uq) 

Quantum mechanical interpretation 

If u(x, t) solves the time-dependent Schrodinger equation 

iut = — u" + x^u 
22 



Harmonic osc. spectrum 
xxxxxx 



Harmonic oscillator modes 

^2 




with potential V(x) = y} and u has XJ norm equal to 1 , then |up represents 
the probability density for the location of a particle in a quadratic potential 
well. 

The k-th eigenfunction Ui^(x) is called the Y-th excited state, because it 
gives a "standing wave" solution 



u(x, t) 



e"'^'^V(x) 



to the time-dependent equation. The higher the frequency or "energy" E^ of 
the excited state, the more it can spread out in the confining potential well, 
as the solution plots show. 

Harmonic oscillator investigations 

Method 1: QDEs Since iio(x) = e^" ^^ is an eigenfunction, we guess that all 
eigenfunctions decay like e^" ^^. So we try the change of variable u = we^" '^^. 
The eigenfunction equation becomes 



w 



2xw'+ (E-1)w = 0, 



which we recognize as the Hermite equation. Solving by power series, we 
find that the only appropriate solutions have terminating power series: they 
are the Hermite polynomials. (All other solutions grow like e" at infinity, 
violating the boundary condition on u.) 

Method 2: Raising and lowering Define 



h- 



dx 



+ X (raising or creation operator) 



dx 



+ X (lowering or annihilation operator). 
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Write H = "1;^+^^ for the harmonic oscillator operator. Then one computes 
that 

H = h+H- + 1 
= h,"h,+ - 1 

Claim. If u is an eigenfunction with eigenvalue E then h^u is an eigenfunc- 
tion with eigenvalue E ± 2. (In other words, h,+ "raises" the energy, and h,^ 
"lowers" the energy.) 
Proof. 

H(H+u) = (h+h" + 1)(h+u) 
= h+(h,"h.+ + l)u 
= h+(H + 2)u 
= h+(E + 2)u 
= (E + 2)h,+u 

and similarly H(h,^u) = (E — 2)h,^u (exercise). 

The only exception to the Claim is that h^u will not be an eigenfunction 
if In u = 0, which occurs precisely when u = Uq = e^" ^^. Thus the lowering 
operator annihilates the ground state. 



Relation to classical harmonic oscillator 

Consider a classical oscillator with mass vn = 2, spring constant k = 2, and 
displacement x(t), so that 2x = — 2x. The total energy is 

x^ + x^ = const. = E. 



To describe a quantum oscillator, we formally replace the momentum x with 



the "momentum operator" — i4- and let the equation act on a function u: 






u = Eu. 



This is exactly the eigenfunction equation — u" + x^u = Eu. 
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Harmonic oscillator in higher dimensions 



— Au + Ixl u = Eu 



Here |xp = xf + ■ ■ ■ +x^. The operator separates into a sum of 1 dimensional 
operators, and hence has product type eigenfunctions 



U = Uk,(Xi)---Uk_j(Xd), 



(2ki + l) + --- + (2ka+l). 



Hydrogen atom in 3 dimensions 



— Au — Au = Eu 



Here V(x) = — 2/|x| is an attractive electrostatic ("Coulomb") potential cre- 
ated by the proton in the hydrogen nucleus. (Notice the gradient of this 
potential gives the correct |x|^^ inverse square law for electrostatic force.) 

Boundary conditions: u(x) ^ as |x| — > oo (we will say more later about 
the precise formulation of the eigenvalue problem) . 

Eigenvalues: E = — 1,— ^,— ^,... with multiplicities 1 ,4, 9, . . . 
That is, the eigenvalue E = — l/n^ has multiplicity n^. 

Eigenfunctions: e^'''^^L^(r)Y{^(9, cj)) for < |Ta| < n— 1, where Y^^ is a 
spherical harmonic and L^ equals r^ times a Laguerre polynomial. 

(Recall the spherical harmonics are eigenfunctions of the spherical Lapla- 
cian in 3 dimensions, with — AgphercY^ = ^{^ + IJYj™'- In 2 dimensions the 
spherical harmonics have the form Y = cos(k9) and Y = sin(k9), which 



satisfy 

Examples 



d^ Y 



k^Y.l 



The first three purely radial eigenfunctions (£ = ra = 0, n 



1,2,3) are e-^ 6-^/^(1 -^), 6-^/^(1 



3 ' ^ 27 ' J 



Hydrogen atom spectilum 
without multiplicitie s 



-1 



4 9 



Hydrogen atom: radial n=l 



u(x)=e 



-^-M 



V(x)=-2/|x| 



The corner in the graph of the eigenfunction at r 
singularity of the Coulomb potential. 



is caused by the 
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Hydrogen atom: radial n=2 



u(x)=e-l^l/^(l — |x|) 
2 




V(x)=-2/|x| 



Hydrogen atom: radial n=3 



2 2 
u<x)=e-W'-'<l — |x|+— |x|^) 

3 27 



V(x)=-2/|x| 



Continuous spectrum Eigenfunctions with positive energy E > do ex- 
ist, but they oscillate as |x| -^ oo, and thus do not satisfy our boundary 
conditions. They represent "free electrons" that are not bound to the nu- 



cleus. See our later discussion of continuous spectrum, in Chapter 18 
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Chapter 4 

Discrete spectral theorem 



Goal 

To state the spectral theorem for an elliptic sesquilinear form on a dense, 
compactly imbedded Hilbert space, and to prove it using the spectral theorem 
for compact, selfadjoint operators. In later chapters we apply the spectral 
theorem to unify and extend the examples of Chapters |2] and [3j 



References BlanchardBriining| Section 6.3 



Matrix preview — weak eigenvectors 

Consider a Hermitian matrix A, and suppose u is an eigenvector with eigen- 
value y, so that Au = yu. Take the dot product with an arbitrary vector v 
to obtain 

Au ■ V = yu ■ V, Vv G C^. 

We call this condition the "weak form" of the eigenvector equation. Clearly 
it implies the original "strong" form, because if ( Au — yu) ■ v = for all v, 
then Au — yu = 0. 

We will find the weak form useful below, when we generalize to Hilbert 
spaces. 

The Hermitian nature of the matrix guarantees conjugate-symmetry of 
the left side of the weak equation, when u and v are interchanged: 



Au ■ V = Av ■ u. 
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This symmetry ensures that all eigenvalues are real, by choosing v = u. We 
will need a similar symmetry property in the Hilbert space setting. 



PDE preview — weak eigenfunctions 

Consider the eigenfunction equation — Au = Au for the Laplacian, in a do- 
main CI. Multiply by a function v G Hj(ri), so that v equals on 90, and 
integrate to obtain 

r 

— vAu dx = A uv dx. 
Jci Jci 

Assume u and v are real-valued. Green's theorem and the boundary condi- 
tion on V imply 

Vu- Vvdx = A(u,v)L2(Q), VvGHj(a). 
a 

We call this condition the "weak form" of the eigenfunction equation. To 
prove existence of ONBs of such weak eigenfunctions, we will generalize to a 
Hilbert space problem. 

Notice the left side of the weak eigenfunction equation, J^ Vu ■ Vv dx, is 
symmetric with respect to u and v. 



Hypotheses 

Consider two infinite dimensional Hilbert spaces "H and /C over M (or C). 
T-L: inner product (u,v)^, norm ||u||^ 
/C: inner product (u, v)^:, norm ||u||x: 

Assume: 

1. /C is continuously and densely imbedded in "H, meaning there exists a 
continuous linear injection l : /C ^ "H with l(/C) dense in "H. 

2. The imbedding K, "-^ Ti is compact, meaning if B is a bounded subset 
of /C then B is precompact when considered as a subset of "H. (Equivalently, 
every bounded sequence in /C has a subsequence that converges in "H.) 

3. We have a map a : /C x /C ^ M (or C) that is sesquilinear, continuous, 
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and symmetric, meaning 

u 1-^ a(u, v) is linear, for eacli fixed v, 

V 1-^ a(u, v) is linear (or conjugate linear), for each fixed u, 

|a(u,v)| < (const.) ||u||;c||v||;c 



a(v, u) = a(u,v) (or a(u,v)) 
4. a is elliptic on /C, meaning 

a(u,u) > c||u||^ Vu G /C, 

for some c > 0. Hence a(u, u) x ||u|||;. 

An important consequence of symmetry and ellipticity is that: 

a(u,v) defines an inner product whose norm is equivalent to the 
||-||x:-norm. 

Spectral theorem 

Theorem 4.1. Under the hypotheses above, there exist vectors Ui , U2, U3, . . . G 
/C and numbers 

< yi < 72 < 73 < ■ ■ ■ ^ 00 

such that: 

• Uj is an eigenvector of a{-,-) with eigenvalue y^, meaning 

a(uj,v) =7j(uj,v)^ VvG/C, (4.1) 

• {uj} is an ONB for T-L, 

• {uj/y^} is an ONB for K. with respect to the a-inner product. 

The decomposition 

f = ^(f,Uj)^u^ (4.2) 

j 

converges in "H for each f G "H, and converges in /C for each f G /C. 
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The idea is to show that a certain "inverse" operator associated with a 
is compact and self adjoint on l-i. This approach makes sense in terms of 
differential equations, where a would correspond to a differential operator 
such as —A (which is unbounded) and the inverse would correspond to an 
integral operator (— A)^^ (which is bounded, and in fact compact, on suit- 
able domains). Indeed, we will begin by solving the analogue of — Au = f 
weakly in our Hilbert space setting, with the help of the Riesz Representation 
Theorem. 

Remark. For a more general spectral theorem, readers may consult the recent 
paper |Auchnnity] and certain references therein. Briefly, the eigenvectors 
there satisfy a(uj,v) = yjb(uj,v) for all v G /C, where the bilinear form b 
is assumed to be weakly continuous on /C x /C. In our situation, b(u, v) = 
(u, v)^ and our assumption that /C imbeds compactly into 1-L implies weak 
continuity of b. 



Proof of Theorem 4-1 We first claim that for each f G 7/ there exists a 



unique u G /C such that 

a(u,v) = (f,v)^ VvG/C. (4.3) 

Furthermore, the map 

f 1-^ u 

is linear and bounded. To prove this claim, fix f G "H and define a bounded 
linear functional F(v) = (v, f)-^ on /C, noting for the boundedness that 

\Hv)\ < ||v||^||f||^ 

< (const.) ||v||K:||f 11^ since K. is imbedded in T-L 



< (const. )a(v,v)^/^||f| 



n 



by ellipticity. Hence by the Riesz Representation Theorem on K, (with respect 
to the a-inner product and norm on /C), there exists a unique u G /C such 
that F(v) = a(v,u) for all v G /C. That is, 

(v,f)w = a(v,u) VvG/C, 



as desired for (4.3). Thus the map B : f i-^ u is well defined. Clearly it is 
linear. And 



a(u,u) = |F(u)| < (const. )a(u,u)^^^||f| 
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n- 



Hence a(u, u)^/^ < (const.) ||f||w, so that B is bounded from T-i to /C, which 
proves our initial claim. 

Next, B : ?^ — > /C — > 7/ is compact, since /C imbeds compactly into l-i. 
Further, B is selfadjoint on "H, since for all f , g G "H we have 



(Bf,g)^ = (g,Bf)H 



= a(Bg, Bf) by definition of B, 

= a(Bf, Bg) by symmetry of a, 

= (f, Bg)^ by definition of B, 

which implies B* = B. 

Hence the spectral theorem for compact, self-adjoint operators [Evans 
App. D] provides an ONB for T-i consisting of eigenvectors of B, with 

Buj = YjUj 



for some eigenvalues y, -^ 0. The decomposition (4.2) holds in T-L because 
{iLj} forms an ONB for "H. 

The eigenvalues of B are all nonzero, because B is injective: Bf = would 



imply (f, v)^ = for all v G /C by (4.3), so that i — (using density of /C in 
n). Thusfj ^0. 

Since we may divide by the eigenvalue we deduce that Uj = B(uj/yj). 
Thus xtj belongs to the range of B, and so Uj G /C. 

The eigenvalues are all positive, since 

yja(ui,v) = a(Buj,v) = (uj,v)>^ Vv G /C 
and choosing v = Uj G /C and using ellipticity shows that y, > 0. Thus we 

def ~ 

see that the reciprocal numbers < yj = 1 /yj -^ cxd satisfy 

a(uj,v) =yj(uj,v)^ \fv e IC, 



which is (4.1). 



Finally, we have a-orthonormality of the set {uj/y/yj}: 

a(uj,Uk) =yj(uj,Uk)^ 
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This orthonorinal set is complete in /C, because if a(itj,v) = for all j then 
(uj,v)^ = for all j, by ( |4.1[ ), so that v = 0. Therefore each f G /C can be 
decomposed as 

j 



with convergence in /C, and this decomposition reduces to (4.2) because 



Remark. Eigenvectors corresponding to distinct eigenvalues are automat- 
ically orthogonal, since 



(yj-yk)(uj,Uk)>^ =yj(uj,Uk)^-yk(uk,Uj)^ 



a(uj,Uk)-a(uk,Uj) 




by symmetry of a. 



32 



Chapter 5 

Application: ONBs of Laplace 
eigenfunctions 



Goal 

To apply the spectral theorem from the previous chapter to the Dirichlet, 
Robin and Neumann Laplacians, and to the fourth order biLaplacian. 



Laplacian 
Dirichlet Laplacian 

— Au = Au in O 
u = on 90 

O =bounded domain in M'^. 

"H = L^(0), inner product (u,v)]_2 = Ji^uvdx. 

/C = Hj(n) =Sobolev space, which is the completion of C^(n) (smooth 
functions equalling zero on a neighborhood of the boundary) under the inner 
product 



(u,v)h1 



[Vu ■ Vv + uv] dx. 



D 



Density: Hj contains C^, which is dense in L^. 
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Continuous imbedding Hj "-^ L^ is trivial: 

vV2 



|U||]_2 



ib''^ 



< 



r N 1/2 

[|Vup + u2]dx^ 
a 

iJ-IIhi 



Compact imbedding: Hg =^ L^ compactly by the Rellich-Kondrachov 
Theorem |GilbargTrudinge"r , Theorem 7.22]. 
Sesquilinear form: define 



a(u,v) 



r r 

Vu • Vv dx + 

a 



uvdx = (u,v)hi, u,vGHj(n). 



a 



Clearly a is symmetric and continuous on Hj(0). 
Ellipticity: a(u,u) = ||u||^, 



The discrete spectral Theorem 4.1 gives an ONB [vlj] for L^(0) and cor- 
responding eigenvalues which we denote yj = Aj + 1 > satisfying 



(uj,v)hi = (Aj + 1)(uj,v)l2 Vv G HI{C1) 



Equivalently, 



Vitj ■ Vv dx = Aj 



a 



itjvdx VvGHj(a) 



Q 



That is, 



-Auj = Ajitj 



weakly, so that Uj is a weak eigenf unction of the Laplacian with eigenvalue Aj. 



Elliptic regularity theory gives that Uj is C°°-smooth in D. GilbargTrudinger 



Corollary 8.11], and hence satisfies the eigenfunction equation classically. 
The boundary condition iij = is satisfied in the sense of Sobolev spaces 
(since Hj is the closure of C^), and is satisfied classically on any smooth 
portion of 911, again by elliptic regularity. 
The eigenvalues are nonnegative, with 

JnjVUiPdx 
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as we see by choosing v = Uj in the weak formulation. 

Further, A, > because: if Aj =0 then |Vuj| = by the last formula, 



so that Uj = by the Sobolev inequality for Hj GilbargTrudingerl Theorem 



7.10], but Uj cannot vanish identically because it has L^-norm equal to 1. 
Hence 

< Ai < A2 < A3 < ■ ■ 



CXD. 



Aside. The Sobolev inequality we used is easily proved: for u G HUO) 



["u-IIl^ 



u^ dx 



a 



by parts 



= — zx^u- — dx 

< 2(max|x|)||u||i_2||9u/9x^||i_2 



< (const.) ||u||i_2 II Vu||l2 

so that we have a Sobolev inequality 

||u||l2 < (const.) II Vu||i_2 Vu G Hj(n), 

where the constant depends on the domain D.. Incidentally, this Sobolev 
inequality provides another proof that Aj > for the Dirichlet Laplacian. 



Neumann Laplacian 



-Au 
9u 
9u 



(J.U 







in a 



on dD. 



n. =bounded domain in M'^ with Lipschitz boundary. 

n = L^ici] _ 

/C = H^ (O) =Sobolev space, which is the completion of C°°(0) under 
the inner product (u, v)hi (see [GilbargTrudingerl p. 174]). 

Argue as for the Dirichlet Laplacian. The compact imbedding is provided 
by the Rellich-Kondrachov Theorem [GilbargTrudinge"? Theorem 7.26], which 
relies on Lipschitz smoothness of the boundary. 



One writes the eigenvalues in the discrete spectral Theorem 4.1 as yj = 
l^j + 1 > and finds 



Vu; ■ Vv dx 



a 



^^j 



a 



Ujvdx 



VvGH^n), 



(5.1) 
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which imphes that 



weakly (and hence classically). In fact (5.1) says a little more, because it 
holds for all v G H^(0), not just for v G Hq(Q) as needed for a weak 
solution. We will use this additional information in the next chapter to show 
that eigenfunctions automatically satisfy the Neumann boundary condition 
(even though we never imposed it)! 

Choosing v = Uj proves \j.j > 0. The first Neumann eigenvalue is zero: 
\i] = 0, with a constant eigenfunction ui = const, j^ 0. (This constant 
function belongs to H^(Il), although not to Hj(0).) Hence 

= M-i < M-i < H3 < • • • ^ cxD. 

Robin Laplacian 

— Au = pu in n 

9u ^ , 

h cru = on oil 

on 

D. ^bounded domain in M'^ with Lipschitz boundary. 

n = L^ici] 

CT > is the Robin constant. 

The density and compact imbedding conditions are as in the Neumann 
case above. 

Before defining the sesquilinear form, we need to make sense of the bound- 
ary values of u. Sobolev functions do have well defined boundary values. 
More precisely, there is a bounded linear operator (called the trace operator) 
T : H^ [D.) -^ L^idCl) such that 

||Tu||L2(gQ) < T||u||Hr(a) (5-2) 

for some T > 0, and with the property that if u extends to a continuous 
function on 90, then Tu = u on 90. (Thus the trace operator truly captures 
the boundary values of u.) Further, if u G Hj(0) then Tu = 0, meaning that 
functions in Hj "equal zero on the boundary". For these trace results, see 
|Evans[ Section 5.5] for domains with C^ boundary, or EvansGariepy, §4.3] 
for the slightly rougher case of Lipschitz boundary. 

36 



Sesquilinear form: 



a[u,vj = 



Vu ■ Vv dx + cr 



n 



uvdS(x) + 



dCl 



uvdx 



n 



(where u and v on the boundary should be interpreted as the trace values 
Tu and Tv). Clearly a is symmetric and continuous on H^ (Q). 
Ellipticity: a(u, u) > ||u||^i, since a > 0. 



One writes the eigenvalues in the discrete spectral Theorem 4.1 as y, = 
Pj + 1 > and finds 



L 



VU; ■ Vv dx + CT 



UjVdS(x) = Pj 



dCL 



a 



itjvdx VvGH^(a), 



which implies that 



-Auj = PjUj 



weakly and hence classically. For the weak solution here we need (by defini- 
tion) only to use trial functions v G Hj(0) (functions equalling zero on the 
boundary). In the next chapter we use the full class v G H^ (II) to show that 
the eigenfunctions satisfy the Robin boundary condition. 
Choosing v = Uj proves 

_ J^|Vujpdx+gJ3^ufdS(x) ^^ 

using again that ct > 0. Further, Pj > because: if pj = then |Vtxj| = so 
that Uj = const., and this constant must equal zero because Jg^ u? dS(x) = 0; 
but Uj cannot vanish identically because it has L^-norm equal to 1 . Hence 
when a > we have 



< Pl < p2 < P3 < 



oo. 



Negative Robin constant: a < 0. Ellipticity more difficult to prove 
when a < 0. We start by controlling the boundary values in terms of the 
gradient and L^ norm. We have 



r r 

u^ dSfx] < (const. 



dCl 



|Vu||u| dx + (const. 



Q 



u^ dx, 



a 
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as one sees by inspecting the proof of the trace theorem ( |Evanst §5.5] or 
|EvansGariepy[ §4.3]). An apphcation of Cauchy-with-e gives 



u^dS(x) < £||Vu||^2 + C||u||^2 

dCL 

for some constant C = C(£) > (independent of u). Let us choose £ = 

l/2|cr|, so that 

1 , 2 

a(u,u) > :r||u||^, — C|(t| ||u||[_2. 

Hence the new sesquihnear form a(u, v) — a(u, v) + C|cr|(u, v)i_2 is elhptic. 
We apply the discrete spectral theorem to this new form, and then obtain the 
eigenvalues of a by subtracting C|ct| (with the same ONB of eigenfunctions). 

Eigenfunction expansions in the L^ and H^ norms 

The L^-ONB of eigenfunctions {uj} of the Laplacian gives the decomposition 

f = ^(f,Uj)L2Uj (5.3) 

j 

with convergence in the L^ and H^ norms, for all f in the following spaces: 

'Hj(a) for Dirichlet, 
f G s H^ [D.) for Neumann, 
H^ fa) for Robin. 



These claims follow immediately from the discrete spectral Theorem |4.1 
view of our applications above. 



m 



Invariance of eigenvalues under translation, rotation and reflection, 
and scaling under dilation 

Eigenvalues of the Laplacian remain invariant when the domain II is trans- 
lated, rotated or reflected, as one sees by a straightforward change of variable 
in either the classical or weak formulation of the eigenvalue problem. Physi- 
cally, this invariance simply means that a vibrating membrane is unaware of 
any coordinate system we impose upon it. 

Dilations do change the eigenvalues, of course, by a simple rescaling re- 
lation: Aj(tO) = t^^Aj(O) for each j and all t > 0, and similarly for the 
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Neumann eigenvalues. (We can understand this scale factor t^^ physically, 
by recalling that large drums vibrate with low tones.) The Robin eigenval- 
ues rescale the same way under dilation, provided the Robin parameter is 
rescaled to a/t on the domain tD.. 

BiLaplacian — vibrating plates 

The fourth order wave equation (t)tt = — AAcj) describes the transverse vibra- 
tions of a rigid plate. (In one dimension, this equation simplifies to the beam 
equation: 4)tt = — (|)""). After separating out the time variable, one arrives 
at the eigenvalue problem for the biLaplacian: 

AAu = Au in 11. 

We will prove existence of an orthonormal basis of eigenfunctions. For sim- 
plicity, we treat only the Dirichlet case, which has boundary conditions 

u = |Vu| = on dCl. 

(The Neumann "natural" boundary conditions are rather complicated, for 
the biLaplacian.) 

O =bounded domain in M'^ 

/C = HofO) = completion of C^(0) under the inner product 
P d d 

(U, v)h2 = [ ^ ^n^XnVx^Xn + ^ ^n.^X,^ + ^^^ ^X. 

•^^ m,n=l m=l 

Density: Hq contains C^, which is dense in L^. 

Compact imbedding: H^ -^ Hj -^ L^ and the second imbedding is com- 
pact. 

Sesquilinear form: define 



a[u,vj 



[ Y- ^''mXnVx^xn + •Liv] dx, u, V G H^(n). 



^ m,n=l 



Clearly a is symmetric and continuous on Ho(0). 
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Ellipticity: ||u||^2 < (d.+ l)a(u, u), because integration by parts gives 



Q 



Uv„,x.„udx 



m=l ra=l 



m=l 

< afu.uld. 



The discrete spectral Theorem 4.1 gives an ONB {uj} for L (O) and cor- 



responding eigenvalues which we denote yj = Aj + 1 > satisfying 

a(uj,v) = (Aj + 1)(uj,v)l2 Vv g H^(a). 
Equivalently, 



o 



L 

Ta,Ti=l 



;^j)x„.XnVx,,XndX = Aj 



a 



itjvdx VvGH^(a). 



That is, 



^i^) 



'^m'^m'^n'^'n 



AjU^ 



Ta,n=l 



weakly, which says 



AAui 



AjUj 



weakly. Hence Uj is a weak eigenfunction of the biLaplacian with eigenvalue 
Aj. Elliptic regularity gives that itj is C°°-smooth, and hence satisfies the 
eigenfunction equation classically. The boundary condition Uj — |Vxtj| = is 
satisfied in the sense of Sobolev spaces (since Uj and each partial derivative 
(uj)x^ belong to Hq), and the boundary condition is satisfied classically on 
any smooth portion of 90, again by elliptic regularity. 
The eigenvalues are nonnegative, with 

A- = ^-^ — > 

' J^ufdx 

as we see by choosing v = Uj in the weak formulation and writing D^u = 



u. 



■mXnJTa,Tl=l 



for the Hessian matrix. 
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Further, Aj > because: if Aj = then (itj)x^xn — by the last formula, 
so that (uj)x„^ = by the Sobolev inequality for Hj applied to (uj)x,^, and 
hence Uj = by the same Sobolev inequality, which gives a contradiction. 
Hence 

< Ai < A2 < A3 < ■ • ■ ^ oo. 

Compact resolvents 



The essence of the proof of the discrete spectral Theorem AA_ is to show 
that the inverse operator B is compact, which means for our differential 
operators that the inverse is a compact integral operator. For example, in 
the Neumann Laplacian application we see that (— A+ 1)^^ is compact from 
L^(0) to H^(n). So is (— A+ oc)^^ for any positive a, but (x — does not 
give an invertible operator because the Neumann Laplacian has nontrivial 
kernel, with — A(c) = for every constant c. 

Thus for the Neumann Laplacian, the resolvent operator 

R, = (-A-A)-^ 

is compact whenever A is negative. 
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Chapter 6 

Natural boundary conditions 



Goal 

To understand how the Neumann and Robin boundary conditions arise "nat- 
urally" from the weak eigenfunction equation. 

Dirichlet boundary conditions 

are imposed directly by our choice of function space Hq(0), since each func- 
tion in that space is a limit of functions with compact support in Q. 

Neumann boundary conditions 

The weak form of the Neumann eigenequation for the Laplacian, from 
Chapter [5} is: 



Vu • Vv dx= \x 



a 



uvdx VveH^(a). 



(6.1: 



Q 



From this formula we showed that — Au = \x\x weakly and hence classically, 
by using only functions v that vanish on the boundary, meaning v G Hj(Q). 
To deduce the Neumann boundary condition 9u/9n = 0, we will take v 
not to vanish on the boundary. Assume for simplicity that the boundary is 
smooth, so that u extends smoothly to O. Green's formula (integration by 



parts) applied to (6.1) implies that 



-Au)v dx + 



a 



dCl 



9u 
dn 



vdS 



(^u)vdx VvGC~(a). 



a 
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Since — Au = ]iu, we deduce 

9u 



3D 



dn 



vdS = VvGC°°(a). 



One may choose v G C°° [D.] to equal the normal derivative of u on the bound- 

(I 
in L^(90); either way one concludes that 



ary (meaning v|aa = f^), or alternatively one may use density of C°°(n)|g^ 



9u 
9rL 



on 9a, 



which is the Neumann boundary condition. 

Note. If the boundary is only piecewise smooth, then one merely applies 
the above reasoning on the smooth portions of the boundary, to show the 
Neumann condition holds there. 

Robin boundary conditions 

Integrating by parts in the Robin eigenfunction equation 



Vu ■ Vv dx + a 



n 



uv dS = p 



an 



uv dx Vv G H^ (O) 



n 



(that is, applying Green's formula to this equation) and then using that 
— Au = pu gives that 

r 9u 

(— + au)vdS = VvGC°°(a). 

Like above, we obtain the Robin boundary condition 

9u 



dn 



+ au = on 90, 



at least on smooth portions of the boundary. 

BiLaplacian — natural boundary conditions 

Natural boundary conditions for the biLaplacian can be derived similarly 
jChasmanl §5]. They are much more complicated than for the Laplacian. 
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Chapter 7 

Application: ONB of 
eigenfunctions for the 
Laplacian with magnetic field 

Goals 

To apply the spectral theorem from Chapter |4] to the magnetic Laplacian 
(the Schrodinger operator for a particle in the presence of a classical magnetic 
field). 

Magnetic Laplacian 

Take a bounded domain D. in M'^, with d = 2 or d = 3. We seek an ONB 
of eigenfunctions and eigenvalues for the magnetic Laplacian 

(iV + A)^u=(3u inn, 
u = on do., 

where u(x) is complex-valued and 

A : Rd ^ ^d 

is a given bounded vector field. 

Physically, A represents the vector potential, whose curl equals the mag- 
netic field: V x A = B. Note that in 2 dimensions, one extends A = 
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(Ai,A2) to a 3-vector (Ai,A2,0) before taking the curl, so that the field 
B = (0, 0, -g^ — g^) cuts vertically through the plane of the domain. For a 
brief explanation of how the magnetic Laplacian arises from the correspon- 
dence between classical energy functions and quantum mechanical Hamilto- 
nians, see |ReedSimon2t p. 173]. 

Now we choose the Hilbert spaces and sesquilinear form. Consider only 
the Dirichlet boundary condition, for simplicity: 



1-i = L^(0; C) (complex valued functions), with inner product 



(u,v)l2 



uvdx. 



Cl 



fC — HJ(0;C) with inner product 



(u,v)h1 



[Vu- Vv + uv]dx 



a 



Density: /C contains C^, which is dense in \J-. 

Continuous imbedding Hj ^^ \? is trivial, since ||u||]_2 < ||u||hi, and the 



imbedding is compact by the Rellich-Kondrachov Theorem [GilbargTrudinger 
Theorem 7.22]. 

Sesquilinear form: define 



a[u,vj 



(iV + A)u ■ (iV + A)v dx + C 



D 



uvdx, u,v G Ho(n; C), 



o 



with constant C = || AH^oo + \. Clearly a is symmetric and continuous on Hj. 
EUipticity: 



> 



[ |Vu|2 + 2 Re(iVu ■ Au) + |A| Vl^ + C|up ] dx 
n 

[ I Vu|^ - 2| Vu||A||u| + 2|Ap|up + \\vi\^ ] dx 



> 



<■ 1 



1 



klVu|2+:^|u|^]dx 



1 



u 



Hi 
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The discrete spectral Theorem 4.1 gives an ONB {uj} for L^(Q; C) and 



corresponding eigenvalues which we denote yj = |3j + C > satisfying 



Q 



(iV + A)uj ■ (iV + A)v) dx = (3j itjvdx VvGHj(a;C). 
n. 

In particular, 

weakly (and hence classically, assuming smoothness of the vector potential 
A), so that Uj is an eigenfunction of the magnetic Laplacian (iV + A)^ with 
eigenvalue (3j. We have 



The eigenvalues satisfy 



|3i < |32 < (33 < >oo. 



n _ Ja I (-LV + AjUjpdx 

Pi — " 



J^lltjPdx 

as we see by choosing v = Uj in the weak formulation. Hence the eigenvalues 
are all nonnegative. 

In fact |3i > if the magnetic field vanishes nowhere, as we will show by 
proving the contrapositive. If |3i =0 then (iV + A)ui = 0, which implies 
A = — iV log Ui wherever ui is nonzero. Then V x A = wherever ui is 
nonzero, since the curl of a gradient vanishes identically. (Here we assume 
Ui is twice continuously differentiable.) Thus the magnetic field vanishes 
somewhere, as we wanted to show. 

Aside. The preceding argument works regardless of the boundary condi- 
tion. In the case of Dirichlet boundary conditions, one need not assume the 
magnetic field is nonvanishing, because the above argument and the reality 
of A together imply that if |3i = then |ui| is constant, which is impossible 
since Ui = on the boundary. 

Gauge invariance 

Many different vector potentials can generate the same magnetic field. 
For example, in 2 dimensions the potentials 

A=(0,xi), A=(-X2,0), A = -(-X2,Xi), 
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all generate the same (constant) magnetic field: V x A = (0,0, 1). Indeed, 
adding any gradient vector Vf to the potential leaves the magnetic field 
unchanged, since the curl of a gradient equals zero. This phenomenon goes 
by the name of gauge invariance. 

How is the spectral theory of the magnetic Laplacian affected by gauge 
invariance? The sesquilinear form definitely changes when we replace A with 
A + Vf . Fortunately, the new eigenf unctions are related to the old by a 
unitary transformation, as follows. Suppose f is C^ -smooth on the closure of 
the domain. For any trial function u G Hj(0; C) we note that the modulated 
function e^^u also belongs to Hj(n;C), and that 



(iV + A)u = (iV + A + Vf ) (e^'u). 

Thus if we write a for the original sesquilinear form and a for the analogous 
form coming from the vector potential A + Vf , we deduce 

a(u,v) = a(e^*^u, e%) 

for all trial functions u,v. Since also (u,v)i_2 = (e^^u, e^^v)i_2, we find that 
the ONB of eigenfunctions Uj associated with a transforms to an ONB of 
eigenfunctions e'-^Uj associated with a. The eigenvalues (energy levels) (3j 
are unchanged by this transformation. 

For geometric invariance of the spectrum with respect to rotations, re- 
fiections and translations, and for a discussion of the Neumann and Robin 



situations, see [LaugesenLiangRoy] Appendix A]. 



Higher dimensions 

In dimensions d > 4 we identify the vector potential A : M'^ ^ M'^ with a 
1 -form 

A = Ai dxi H h Aa dx^ 

and obtain the magnetic field from the exterior derivative: 

B = dA. 
Apart from that, the spectral theory proceeds as in dimensions 2 and 3. 
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Chapter 8 

Application: ONB of 
eigenfunctions for Schrodinger 
in a confining well 



Goal 

To apply the spectral theorem from Chapter |4] to the harmonic oscillator and 
more general confining potentials in higher dimensions. 

Schrodinger operator with potential growing to infinity 

We treat a locally bounded, real-valued potential V(x) on M'^ that grows 
at infinity: 

— C < V(x) ^ oo as |x! -^ cxD, 

for some constant C > 0. For example, V(x) = |xp gives the harmonic 
oscillator. 

We aim to prove existence of an ONB of eigenfunctions and eigenvalues 
for 

[-A + V)u = Eu in R'^ 

u ^ as Ixl ^^ oo 



V. = L^(]R'^), inner product (u,v)l2 = Jj^a uv dx. 
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/C = H^ [R'^) n L^(|V| dx) under the inner product 



(u,v)^ 



[Vu-Vv+(1 +|V|)uv]dx. 



Density: /C contains C^, which is dense in L^. 

Continuous imbedding /C '^^ L^ is trivial, since ||u||i_2 < ||u||x:. To prove 
the imbedding is compact: 

Proof that imbedding is compact. Suppose {fi<} is a bounded sequence in /C, 
say with llfklk < M. for all k. We must prove the existence of a subsequence 
converging in L^(]R'^). 

The sequence is bounded in H^ (B(R)) for each ball B(R) C M'^ that is cen- 
tered at the origin. Take R = 1 . The Rellich-Kondrachov theorem provides 
a subsequence that converges in L^(B(1)). Repeating with R = 2 provides a 
sub-subsequence converging in L^(B(2)). Continue in this fashion and then 
consider the diagonal subsequence, to obtain a subsequence that converges 
in L2(B(R)) for each R > 0. 

We will show this subsequence converges in L^(]R'^). Denote it by {fkj- 
Let £ > 0. Since V(x) grows to infinity, we may choose R so large that 
V(x) > 1/e when |x| > R. Then 



Rd\B(R) 



fl dx < £ 



Kd\B(R) 



f'.Vdx 



< pllf l|2 
S £||Tkelk 



for all L Since also {fkJ converges on B(R), we have 
fknJk^fKdi = hmsupllfv,, — f- 



limsup||fi^j 



limsupllf,,^ - fk„J|L2(Rd\B(R)) < 2^/eM. 



Therefore {fkJ is Cauchy in L^lM^), and hence converges. 
Sesquilinear form: define 



[Vu- Vv + Vuv]dx+ (2C + 1) uvdx, u,vg/C 



Clearly a is symmetric and continuous on /C. 

Ellipticity: a(u,u) > ||u||^, since V + 2C+ 1 > 1 + |V|. 
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The discrete spectral Theorem 4.1 gives an ONB {itj} for L^flR'^) and 
corresponding eigenvalues which we denote yj = Ej + 2C + 1 > satisfying 



[ Vxtj ■ Vv + Vujv] dx 



UjV dx Vv G /C. 



EjUj 



In particular, 

-Auj + Vuj 

weakly (and hence classically, assuming smoothness of V), so that Uj is an 
eigenfunction of the Schrodinger operator —A + V, with eigenvalue Ej. We 
have 

ti < E2 < E3 < > 00. 

The boundary condition Uj ^ at infinity is interpreted to mean, more 
precisely, that u, belongs to the space H^(]R'^) fl L^(|V| dx). This condition 
suffices to rule out the existence of any other eigenvalues for the harmonic 
oscillator, for example, as one can show by direct estimation [Strauss] . 

The eigenvalues satisfy 

' JlRd Uf dx 

as we see by choosing v = Uj in the weak formulation. Hence if V > then 
the eigenvalues are all positive. 
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Chapter 9 

Variational characterizations of 
eigenvalues 



Goal 

To obtain minimax and maximin characterizations of the eigenvalues of the 
sesquilinear form in Chapter |4j 

References |Bandle] Section 111.1.2 

Motivation and hypotheses. How can one estimate the eigenvalues if the 
spectrum cannot be computed explicitly? We will develop two complemen- 
tary variational characterizations of eigenvalues. The intuition for these char- 
acterizations comes from the special case of eigenvalues of a Hermitian (or 
real symmetric) matrix A, for which the sesquilinear form is a(u,v) = Au-v 
and the first eigenvalue is 

Av- V 
yi = mm ^. 

v^O V ■ V 

We will work under the assumptions of the discrete spectral theorem in 
Chapter |4| for the sesquilinear form a. Recall the ordering 

Ti < 72 < 73 < ■ ■ ■ ^ oo. 
Poincare's minimax characterization of the eigenvalues 
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Define the Rayleigh quotient of u to be 

a(u,u) 

We claim yi equals tlie minimuni value of tlie Rayleigli quotient: 



. a(f,f) 
yi = mm 



f6/C\{0} (f , f ) 



H 



(9.i: 



This characterization of the first eigenvalue is the Rayleigh principle. 

More generally, each eigenvalue is given by a minimax formula known as 
the Poincare principle: 



y; = mm max ^ — 

s fes\{0} (f , f )^ 



(9.2) 



where S ranges over all j-dimensional subspaces of /C. 

Remark. The Rayleigh and Poincare principles provide upper bounds on 
eigenvalues, since they expresses y, as a minimum. More precisely, we obtain 
an upper bound on yj by choosing S to be any j-dimensional subspace and 
evaluating the maximum of the Rayleigh quotient over f G S. 

Proof of Poincare principle. First we prove the Rayleigh principle for the 
first eigenvalue. Let f G /C. Then f can be expanded in terms of the ONB of 
eigenvectors as 

f = ^ Cj-L^ 

j 

where Cj = (f, Uj)^. This series converges in both "H and /C (as we proved in 
Chapter |4]). Hence we may substitute it into the Rayleigh quotient to obtain 

a(f,f) Li,kCjC^a(uj,Uk) 



[9.3] 



LjICjPTj 



2_i IHI 



since the eigenvectors {uj} are orthonormal in H and the collection {i tj/^/ Yj} 
is a-orthonormal in /C (that is, a(uj,Uk) = yjSjk)- The expression ( |9.3[ ) is 
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obviously greater than or equal to yi , with equality when f = Ui , and so we 



have proved the Rayleigh principle (9.1). 



Next we prove the minimax formula (9.2) for j = 2. (We leave the case 
of higher j-values as an exercise.) Choose S — {ciUi + C2U2 : Ci , C2 scalars} to 
be the span of the first two eigenvectors. Then 

a(f,f) Li=ilCjl'yj 
max ^ — = max — -^ = 72- 

fes\{0} (f , f )^ (ci,c2)^{o,o) ^^^1 |Cjp 



Hence the minimum on the right side of (9.2[) is < y2. 



To prove the opposite inequality, consider an arbitrary 2-dimensional sub- 
space S C /C. This subspace contains a nonzero vector g that is orthogonal to 
ui (since given a basis {vi,V2} for the subspace, we can find scalars di, d2 not 
both zero such that g = diVi + d2V2 satisfies = di (vi , ui )^ + d2(v2, ui )^ = 



(g,ui)^). Then Ci = in the expansion for g, and so by (9.3), 

a(g,g) ^ lL]=i\^]\i) 
(g>g)w 2_j=2lcjr 

Hence 

a(f,f) ^ a(g,g) 

fes\{0} f,f^ g,gw 



which implies that the minimum on the right side of (9.2) is > yi. D 



Variational characterization of eigenvalue sums. The sum of the 

first n eigenvalues has a simple "minimum" characterization, similar to the 
Rayleigh principle for the first eigenvalue, but now involving pairwise orthog- 
onal trial functions: 

yi + --- + yn (9.4) 

= min|^%lll + ... + ^^:fjG/C\{0},(f^,f,)«=Owhenj^k|. 
See Handle's book for the proof and related results [Bandlet Section III. 1.2]. 
Courant's maximin characterization 
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The eigenvalues are given also by a maximin formula known as the Courant 
principle: 



Yj = max mm 



a(f,f) 



s fes^\{0} (f , f ) 



v. 



(9.5) 



where this time S ranges over all (j — 1 ) -dimensional subspaces of /C. 



Remark. The Courant principle provide lower hounds on eigenvalues, since 
it expresses y, as a maximum. The lower bounds are difficult to compute, 
however, because S"*" is an infinite dimensional space. 

Sketch of proof of Courant principle. The Courant principle reduces to 
Rayleigh's principle when j = 1 , since in that case S is the zero subspace 
and S^ = /C. 

Now take j = 2 (we leave the higher values of j as an exercise). For 
the "<" direction of the proof, we choose S to be the 1 -dimensional space 
spanned by the first eigenvector Ui. Then every f G S"*" has Ci = (f, Ui)^ = 
and so 

. . a(f>f) 
72 < mm 

fesi\{0}(f,f)« 

by expanding f = YJ^i '-j^j ^^'^ computing as in our proof of the Poincare 
principle. 

For the ">" direction of the proof, consider an arbitrary 1 -dimensional 
subspace S of /C. Then S"*" contains some vector of the form f = CiUi + C2U2 
with at least one of Ci or C2 nonzero. Hence 

. a(Vf) L^ikjpTj ^ 
fes^\{0} (f , f )^ Lj=i ICjP 

as desired. 



Eigenvalues as critical values of the Rayleigh quotient 

Even if we did not know the existence of an ONB of eigenvectors we could 
still prove the Rayleigh principle, by the following direct approach. Define 
y* to equal the infimum of the Rayleigh quotient: 

y = mt ^ — . 

fe^\{0} (f , f )^ 
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We will prove y* is an eigenvalue. It follows that y* is the lowest eigenvalue, 
y* = yi (because if any eigenvector f corresponded to a smaller eigenvalue, 
then the Rayleigh quotient of f would be smaller than y*, a contradiction). 
First, choose an infimizing sequence {fi^} normalized with ||fk||w = 1, so 
that 

a(fk,fk) ^y*. 

By weak compactness of closed balls in the Hilbert space /C, we may suppose 
fk converges weakly in /C to some u G /C. Hence fk also converges weakly in 
T-L to u (because if ¥{■) is any bounded linear functional on "H then it is also 
a bounded linear functional on /C). We may further suppose fk converges in 
T-L to some v G "H (by compactness of the imbedding /C "-^ Ti) and then fk 
converges weakly in "H to v, which forces v = u. To summarize: f k ^ u 
weakly in /C and fk — ^ u in "H. In particular, ||u||^ = 1. Therefore we have 

< a(fk — u, fk — u) 
= a(fk, fk) — 2 Re a(fk, u) + a(u, u) 

-^ y* — 2 Re a(u, u) + a(u, u) using weak convergence f k ^ u 
= y*-a(u,u) 
<0 

by definition of y* as an infimum. 

We have shown that the infimum defining y* is actually a minimum, 

. a(f,f) 
y = mm ^ — , 

feAlo} (f , f)n 

and that the minimum is attained when f = u. 

Our second task is to show u is an eigenvector with eigenvalue y*. Let 
V G /C be arbitrary and use f = u + £V as a trial function in the Rayleigh 
quotient; since u gives the minimizer, the derivative at £ = must equal zero 
by the first derivative test from calculus: 



d a(u+£v, u+£v) 



d£ (U+£V,U+£V)^ 



= 2Re a(u,v) — y*2Re(u,v)^. 

£=0 



The same equation holds with Im instead of Re, as we see by replacing v with 
iv. (This last step is unnecessary when working with real Hilbert spaces, of 
course.) Hence 

a(u, v) = y* (u, v)n Vv G /C, 
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which means u is an eigenvector for the sesquihnear form a, with eigenvalue 

y*- 

Aside. The higher eigenvalues (yj for j > 1 ) can be obtained by a similar 
process, minimizing the Rayleigh quotient on the orthogonal complement of 
the span of the preceding eigenf unctions Ui , . . . , xtj_i . In particular, 

. a(f,f) 
yi = mm 



f±ui (f,f) 



H 



where f _L Ui means that the trial function f G /C\{0} is assumed orthogonal 
to Ui in H: (f, Ui)^ = 0. 
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Chapter 10 

Monotonicity properties of 
eigenvalues 

Goal 

To apply Poincare's ininiinax principle to the Laplacian and related opera- 
tors, and hence to establish monotonicity results for Dirichlet and Neumann 
eigenvalues of the Laplacian, and a diamagnetic comparison for the magnetic 
Laplacian. 

References |Bandle] 

Laplacian, biLaplacian, and Schrodinger operators 



Applying the Rayleigh principle (9.1) to the examples in Chapters 5148 gives 



Jo iVfP dx 

A] = min -^ — z Dirichlet Laplacian on Q, 

feHj(a) J^f^dx 

LlVfpdx+aL^f^dS 

pi = min -^ ^ — ^^ Robin Laplacian on CI, 

feH' (Q) J ^ f^ dx 

Jo iVf P dx 
(Xi = min — 7; — z Neumann Laplacian on O, 

f6H' (CI) J ^ f-^ dx 
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Ai = mill — — ^'"'~ ^""'^'^ Dirichlet biLaplacian on Q. 



mm 



Ja(Af)^dx 

Ld li-Vf + Af P dx 
(3i = mill — !^ — 7; — — magnetic Laplacian 

feHi(a) JjgJfKdx 

Ld(lVfp + Vf2)dx 
El = min ^-^ Schrodinger with potential 

feHl (Kd)nL2(|V| dx) J jja f Cix 

V(x) growing to infinity. 

The Poincare principle applies too, giving formulas for the higher eigenvalues 
and hence implying certain monotonicity relations, as follows. 

Neumann < Robin < Dirichlet 

Free membranes give lower tones than partially free and fixed membranes: 

Theorem 10.1 (Neumann-Rob in-Dirichlet comparison). Let CI be a bounded 
domain in M'^ with Lipschitz boundary, and fix a > 0. 

Then the Neumann eigenvalues of the Laplacian lie below their Robin 
counterparts, which in turn lie below the Dirichlet eigenvalues: 

^j<Pj<Aj Vj>l. 

Proof. Poincare's minimax principle gives the formulas 

Jg iVf 1^ dx 
LI; = mm max -^ — 

' u feu\{0} J ^ f 2 dx 

J^lVfpdx+gJa^f^dS 

Pi = mm max ^^-^ -j. — - — ^-^^ 

T feT\{0} J ^ f^ dx 

. L iVf I' dx 

A; = mm max — 7. — 

s fes\{0} J ^ f^ dx 

where S ranges over all j-dimensional subspaces of Hj(ri), and T and U range 
over all j-dimensional subspaces of H^ (O). 

Clearly fJ^ < Pj. Further, every subspace S is also a valid T, since Hj C H^ 
Thus the minimum for Pj is taken over a larger class of subspaces. Since for 
f G HJ the boundary term vanishes in the Rayleigh quotient for Pj, we 
conclude that Pj < Aj. D 
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Domain monotonicity for Dirichlet spectrum 

Making a drum smaller increases its frequencies of vibration: 

Theorem 10.2. Let D. and D. be bounded domains in M"^, and denote the 
eigenvalues of the Dirichlet Laplacian on these domains by Aj and Aj, respec- 
tively. 

IfCLDD. then 

Aj < Aj Vj > 1 . 

Proof. Poincare's minimax principle gives that 

. Jg iVf P dx 
Aj = mm max -^ — 

' s fes\{0} J ^ f^ dx 

T . I a iVf P dx 

Ai = mm max 



s fes\{0} 



Jof'dx 



where S ranges over all j-dimensional subspaces of Hj(Il) and S ranges over 
all j-dimensional subspaces of Hj(Il). 

Every subspace S is also a valid S, since Hj(0) C Hj(n) (noting that 
any approximating function in C^(0) belongs also to 0^(0) by extension 
by 0.) Therefore Aj < Aj. D 

Restricted reverse monotonicity for Neumann spectrum 

The monotonicity proof breaks down in the Neumann case because H^ [Q.) is 
not a subspace of H^ (O). More precisely, while one can extend a function in 
H^ (O) to belong to H^ (Q), the extended function must generally be nonzero 
outside n, and so its L^ norm and Dirichlet integral will differ from those of 
the original function. 

Furthermore, counterexamples to domain monotonicity are easy to con- 
struct for Neumann eigenvalues, as the figure below shows with a rectangle 
contained in a square. In that example, the square has side length 1 and hence 
[iz = 7t^, while the rectangle has side length \/2[0.9) and so JI2 = 7r^/(1.62), 
which is smaller than (j.2- 

Nonetheless, monotonicity does holds in a certain restricted situation, 
although the inequality is reversed — the smaller drum has lower tones: 

59 




Theorem 10.3. Let O and O. be bounded Lipschitz domains in M.^, and 
denote the eigenvalues of the Neumann Laplacian on these domains by (x,- 
and ]Ij, respectively. 

If D. G D. and Q.\D. has measure zero, then 

Hj < ^j Vj > 1 . 

One might imagine the smaller domain O as being constructed by remov- 
ing a hypersurface of measure zero from Q, thus introducing an additional 
boundary surface. Reverse monotonicity then makes perfect sense, because 
the additional boundary, on which values are not specified for the eigenfunc- 
tions, enables the eigenf unctions to "relax" and hence lowers the eigenvalues. 





Introducing additional boundary surfaces to a Dirichlet problem would 
have the opposite effect: the eigenfunctions would be further constrained, 
and the eigenvalues raised. 
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Proof. Poincare's ininimax principle gives that 

s"' fes\{0} J^ f ^ dx 

L I vf I' dx 

s fes\{0} ici^^'^^ 



Lij = mm max 

' s fes\{o: 

]Ij = mill miax 



where S ranges over all j- dimensional subspaces of H^ [D.) and S ranges over 
all j-dimensional subspaces of H^ (Q). 

Every subspace S is also a valid S, since each f G H^ (O) restricts to a 
function in H^ (Q) that has the same H^-norm (using here that D.\Q. has 
measure zero). Therefore \ij < \Xj. D 

Diamagnetic comparison for the magnetic Laplacian 

Imposing a magnetic field always raises the ground state energy. 
Theorem 10.4 (Diamagnetic comparison). 

|3i >Ai 
First we prove a pointwise comparison. 
Lemma 10.5 (Diamagnetic inequality). 

l(iV + A)f|> |V|f|| 
Proof of Lemma\lO.^ Write f in polar form as f = Re^®. Then 



|iVf + Af 1^ = lie^^VR - Re^^Ve + ARe^^l^ 
= |iVR - RVe + AR|^ 

= |VR|^ + R^|VG-A|^ 
> |VR|^ = |V|f||^ 



n 



Proof of Theorem 10. 4 The proof is immediate from the diamagnetic in- 



equality in Lemma |10.5| and the Rayleigh principles for |3i and Ai at the 
beginning of this chapter. Note we can assume f > in the Rayleigh prin- 
ciple for Ai , since the first Dirichlet eigenfunction can be taken nonnegative 



[GilbargTrudingeT Theorem 8.38]. D 



61 



Chapter 11 

Weyl's asymptotic for high 
eigenvalues 

Goal 

To determine the rate of growth of eigenvalues of the Laplacian. 

References jArendt j : |CourantHilbert] Section VI. 4 

Notation 

The asymptotic notation ct, ~ (3j means 

lim — ^ = 1 . 
Write Vd for the volume of the unit ball in d-dimensions. 

Growth of eigenvalues 

The eigenvalues of the Laplacian grow at a rate cj^'^'^ where the constant 
depends only on the volume of the domain, independent of the boundary 
conditions. 
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Theorem 11.1 (Weyl's law). Let D. be a bounded domain in M'^ with piece- 
wise smooth boundary, ^s j ^ oo the eigenvalues grow according to: 



(i7lj/\Cl\]' 

Aj ~ pj ~ ^j ~ <^ 4n]/\a\ 

i(67r^j/|a|)'/' 



;d 
d 



1) 

2) 
3) 



and more generally, 



Aj ~ pj ~ |ij ~ 47t^ 



] 



VdlOl 



2/d 



:d> r 



Here |0| denotes the &- dimensional volume of the domain, in other words its 
length when d — 1 and area when d = 2. 

In 1 dimension the theorem is proved by the exphcit formulas for the 
eigenvalues in Chapter [2] We will prove the theorem in 2 dimensions, by 
a technique known as "Dirichlet-Neumann bracketing". The higher dimen- 
sional proof is similar. 

An alternative proof using small-time heat kernel asymptotics can be 
found (for example) in the survey paper by Arendt et al. [Arendtt §1.6]. 

Proof of Weyl aymptotic — Step 1: rectangular domains. In view of the Neu- 



mann-Robin-Dirichlet comparison (Theorem 10.1 ), we need only prove Weyl's 
law for the Neumann and Dirichlet eigenvalues. We provided a proof in 
Proposition 2.1[ for rectangles. D 



Proof of Weyl aymptotic — Step 2: finite union of rectangles. Next we sup- 
pose Ri , . . . , Rn are disjoint rectangular domains and put 

n = u-^iR^,_ 
a = int {ui^^,^:) . 



n 



n 
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For example, if Ri and R2 are adjacent squares of side length 1 , then Q. is 
the disjoint union of those squares whereas O is the 2 x 1 rectangular domain 
formed from the interior of their union. 

Admittedly D. is not connected, but the spectral theory of the Laplacian 
remains valid on a finite union of disjoint domains: the eigenf unctions are 
simply the eigenfunctions of each of the component domains extended to 
be zero on the other components, and the spectrum equals the union of 
the spectra of the individual components. (On an infinite union of disjoint 
domains, on the other hand, one would lose compactness of the imbedding 
H^ ^^ L^, and the zero eigenvalue of the Neumann Laplacian would have 
infinite multiplicity.) 

Write Aj and ]Ij for the Dirichlet and Neumann eigenvalues of O. 



Then by the restricted reverse Neumann monotonicity (Theorem 10.3) 



Neumann-Robin-Dirichlet comparison (Theorem 10.1) and Dirichlet mono- 



tonicity (Theorem 10.2), we deduce that 



l^j < ^j < Pj < Aj < Aj 
Hence if we can prove Weyl's law 



Mi~Aj 



47tj 



Vj > 1. 



11.1^ 



for the union-of-rect angles domain II, then Weyl's law will follow for the 
original domain D.. 

Define the eigenvalue counting functions of the rectangle Rra to be 

NNeu('x;R^) = #{j > 1 : [l^[R^] < a}, 
NDir(cx;R^) = #{j > 1 : Aj(R^) < a}. 

We know from Weyl's law for rectangles (Step 1 of the proof above) that 



NNeu(a;RTu) ~ NDir(a;Rr 



IR. 



47t 



-(X 



11.21 



as a ^ cxD. 

The spectrum of D. is the union of the spectra of the Rm,, and so (here 
comes the key step in the proof!) the eigenvalue counting functions of O 
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equal the sums of the corresponding counting functions of the rectangles: 

n 

NNcu(a;n) = Y NNcu(oc;RTn), 



m=l 



N 



Dir ( (x; Q) = ^ N Dir ( oc; K 



m=l 



Combining these sums with the asymptotic (|11.2) shows that 

ini 



NNeu(a;n) 




47t 



-a 



and similarly 



NDir(a;n) 



4n 



a 



as (X ^ oo . W e can invert these last two asymptotic^ formulas with the help 
of Lemma 



2.2 



thus obtaining Weyl's law (11.1) for D. 



a 



Proof of Weyl aymptotic — Step 3: approximation of arbitrary domains. Lastly 
we suppose O is an arbitrary domain with piecewise smooth boundary. The 
idea is to approximate D. with a union-of-rectangles domain such as in Step 
2, such that the volume of the approximating domain is within e of the vol- 
ume of II. We refer to the text of Courant and Hilbert for the detailed proof 
|CourantHilb"irtl §VI.4.4]. D 
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Chapter 12 

Polya's conjecture and the 
Berezin— Li— Yau Theorem 

Goal 

To describe Polya's conjecture about Weyl's law, and to state the "tiling 
domain" and "summed" versions that are known to hold. 



References [Kellnert [Laptev P61ya| 



Polya's conjecture 



Weyl's law (Theorem 11.1) says that 



, 47tj 

^ ~ M ~ ^' as 3 ^ oo, 

for a bounded plane domain O with piecewise smooth boundary. (We restrict 
to plane domains, in this chapter, for simplicity.) 

Polya conjectured that these asymptotic formulas hold as inequalities. 

Conjecture 12.1 ( |P61ya| , 1960). 

4™ 
Ai>^>l^j Vj>l. 

The conjecture remains open even for a disk. 

Polya proved the Dirichlet part of the inequality for tiling domains |P61ya| , 
and Kellner did the same for the Neumann part [Kellnerj . Recall that a 
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"tiling domain" covers the plane with congruent copies of itself (transla- 
tions, rotations and reflections). For example, parallelograms and triangles 
are tiling domains, as are many variants of these domains (a fact that M. C 
Escher exploited in his artistic creations). 

Polya and Kellner's proofs are remarkably simple, using a rescaling argu- 
ment together with Weyl's law. 

For arbitrary domains, Polya's conjecture has been proved only for Ai, Ai 
(see [Henrott Th. 3.2.1 and (4.3)]) and for |a.i, |j.2, 1^3 (see [Girouard] ) . The 
conjecture remains open for j > 3 (Dirichlet) and j > 4 (Neumann). 

Berezin— Li— Yau results 

The major progress for arbitrary domains has been on a "summed" version of 
the conjecture. (Quite often in analysis, summing or integrating an expres- 
sion produces a significantly more tractable quantity.) Li and Yau |LiYau] 
proved that 

2nf 



z 



Ak> 



which is only slightly smaller than the quantity (27t/|£l|)j(j + 1 ) that one gets 
by summing the left side of the Polya conjecture. An immediate consequence 
is a Weyl-type inequality for Dirichlet eigenvalues: 

by combining the very rough estimate jAj > Y-k=^ ^^ with the Li- Yau in- 
equality. The last formula has 2n whereas Polya's conjecture demands An, 
and so we see the conjecture is true up to a factor of 2, at worst. 

Similar results hold for Neumann eigenvalues. 

A somewhat more general approach had been obtained earlier by Berezin. 
For more information, consult the work of Laptev [Laptev] and a list of open 
problems from recent conferences |AIMj . 
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Chapter 13 

Case study: stability of steady 
states for reaction— diffusion 
PDEs 



Goal 

To linearize a nonlinear reaction-diffusion PDE around a steady state, and 
study the spectral theory of the linearized operator by time-map methods. 

References [Schaaf] Section 4.1 

Reaction— diffusion PDEs 

Assume throughout this section that ffy) is a smooth function on M. Let 
X > 0. We study the reaction-diffusion PDE 

Ut = -a.. + f(u) (13.1) 

on the interval (0, X) with Dirichlet boundary conditions u(0) = u(X) = 0. 
Physical interpretations include: (i) u ^temperature and f =rate of heat 
generation, (ii) u ^chemical concentration and f ^reaction rate of chemical 
creation. 

Intuitively, the 2nd order diffusion term in the PDE is stabilizing (since 
xtt = Uxx is the usual diffusion equation), whereas the 0th order reaction 
term can be destabilizing (since solutions to Ut = f (u) will grow, when f is 
positive). Thus the reaction-diffusion PDE features a competition between 
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stabilizing and destabilizing effects. This competition leads to nonconstant 
steady states, and interesting stability behavior. 

Steady states. If U(x) is a steady state, then 

U" + f (U) =0, < X < X. (13.2) 

More than one steady state can exist. For example if f (0) = then U = 
is a steady state, but nonconstant steady states might exist too, such as 
U(x) = sinx when X = 7t and f(ij) = y. 

Linearized PDE 

We perturb a steady state by considering 

U = U+ £({) 

where the perturbation 4)(x, t) is assumed to satisfy the Dirichlet BC cf) = 



at X = 0, L, for each t. Substituting u into the equation (13.1) gives 

+ £(t)t = (U^x + £4)xx) + f (U + ecf)) 

= U,, + £Ct),, + f (U) + f (U)£4) + 0(£^). 

The leading terms, of order £°, equal zero by the steady state equation for 
U. We discard terms of order £^ and higher. The remaining terms, of order 
e\ give the linearized equation: 

4)t = ct).. + f'(U)4). (13.3) 

That is, 

where L is the symmetric linear operator 

Lw = — Wxx — f'(U)w. 
Separation of variables gives (formally) solutions of the form 

i 
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where the eigenvalues Tj and Dirichlet eigenfunctions Wj satisfy 



LWj = TjWj 

with Wj(0) =Wj(X) =0. 

Thus the steady state U of the reaction-diffusion PDE is 

linearly unstable if Ti < 

because the perturbation cj) grows to infinity, whereas the steady state is 

linearly stable if Ti > 

because ({) remains bounded in that case. 

To make these claims rigorous, we study the spectrum of L. 

Spectrum of L 

We take: 

n=(o,x) 

H = L^(0, X), inner product (u,v)]_2 = J^ uv dx 
/C = Hq(0, X), inner product 



rX 



(u,v)hi 



fuV + uv) dx 



Compact imbedding Hj ^^ \? by Rellich-Kondrachov 
Symmetric sesquilinear form 



(uV-f'(U)uv+Cuv)dx 



where C > is chosen larger than ||f' ||l°° + 1- Proof of ellipticity: 



a[u,uj > 



((u')^ + u^)dx= ||u| 



2 
Hi 



by choice of C. 



The discrete spectral Theorem 4.1 now yields an ONB of eigenfunctions 
{wj} with eigenvalues y, such that 

a(wj,v) =yj(wj,v)i_2 VvgHJ(0,X). 



70 



Writing yj = Tj + C we get 

pX pX 











Wjvdx VvgHJ(0,X). 



(wjV - f'(U)wjv) dx = Tj 
These eigenfunctions satisfy Lwj = TjWj weakly, and hence also classically. 

Stability of the zero steady state. 

Assume f (0) = 0, so that U = is a steady state. Its stabihty is easily 
determined, as follows. 

The linearized operator is Lw = — w" — f{0)w, which on the interval 
(0, X) has Dirichlet eigenvalues 

^i = (|)'-m. 

Thus the zero steady state is linearly unstable if and only if 

Thus we may call the reaction-diffusion PDE "long-wave unstable" when 
f'(0) > 0, because then the zero steady state is unstable with respect to per- 
turbations of sufficiently long wavelength X. On short intervals, the Dirichlet 
BCs are strong enough to stabilize the steady state. 

Sufficient conditions for linearized instability of nonconstant steady 
states 

Our first instability criterion is structural, meaning it depends on properties 
of the reaction function f rather than on properties of the particular steady 
state U. 

Theorem 13.1. Assume the steady state U is nonconstant, and that f (0) = 
0, f'iO] = and f" > 0. Then Ti < 0. 

For example, the theorem shows that nonconstant steady states are un- 
stable when f (y) = y^. 
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Ti = mill i ° '' — '-^ ^ — : w G 



Hj(0,X)}. 



Proof. First we collect facts about boundary values, to be used later in the 
proof when we integrate by parts: 

U = at X = 0, X by the Dirichlet BC, 

f (U) = at X = 0, X since f (0) = 0, 

U" = at X = 0, X because U" = -f (U) 

f" (U) = at X = 0, X since f"{0] = 0. 

The Rayleigh principle for L says that 

Jo((w -)^-f(U) w^)dx 
J^w^dx 

We choose a trial function 

w = U", 

which is not the zero function, since U is nonconstant. Then the numerator 
of the Rayleigh quotient for w is 

pX 

((U'")2_f'(u)(U")')dx 



rX 

= ( - U"" - f'(U)U") U" dx by parts 
Jo 





1 



f"(U) (U')^U" dx by the steady state equation ( |13.2| ) 

rX 



f"(U)[(U')^]'dx 
f'lUJlUO'^dx by parts 



^ Jo 

3 Jo 
<0 

since f" > and U is nonconstant. Hence Ti < 0, by the Rayleigh principle. 

D 

Motivation for the choice of trial function. Our trial function w — U" corre- 
sponds to a perturbation u = U + £4) = U+ ell", which tends (when e > 0) 
to push the steady state towards the constant function. The opposite per- 
turbation (e < 0) would tend to make the solution grow even further away 
from the constant steady state. 
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The next instability criterion, rather than being structural, depends on 
particular properties of the steady state. 

Theorem 13.2 ( jSchaaf I Proposition 4.1.2]). Assume the nonconstant steady 
state U changes sign on (0,X). Then x-\ < 0. 

For example, suppose f (y) = "y so that the steady state equation is U" + 
U = 0. If X = In then the steady state U — sinx is linearly unstable, by 
the theorem. Of course, for that example we can compute the spectrum of 
L exactly: the lowest eigenfunction is w = sin(x/2) with eigenvalue Ti = 

(1)^-1 <0. 

Proof. If U changes sign then it has a positive local maximum and a negative 
local minimum in (0, X), recalling that U = at the endpoints. Obviously 
U' must be nonzero at some point between these local extrema, and so there 
exist points <x-\ < xi < X such that 

U'(xi) =U'(X2) =0 

and U' 7^ on (xi,X2). Define a trial function 

U' on (xi.Xi), 
w = < 

elsewhere. 

(We motivate this choice of trial function at the end of the proof.) Then 
w is piecewise smooth, and is continuous since w = U' = at Xi and Xz- 
Therefore w G Hj(0, X), and w ^ since U' ^ on (xi,X2). 
The numerator of the Rayleigh quotient for w is 

{[w']^-f[U]w^)dx= (-w"-f'(U)w)wdx by parts 



smce 







w 



-U'" = (f(U))' = f (U)U' = f (U)w. (13.4) 



Hence Ti < 0, by the Rayleigh principle for the first eigenvalue. 

Suppose Ti = 0, so that the Rayleigh quotient of w equals Ti. Then 
w must be an eigenfunction with eigenvalue Ti (because substituting w = 
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Y_; CjWj into the Rayleigh quotient would give a value larger than Ti, if Cj 
were nonzero for any term with eigenvalue larger than Ti). 

Since the eigenfunction w must be smooth, the slopes of w from the left 
and the right at X2 must agree, which means yv'lxi) = w'(x2+) = 0. Thus 



"^^(^2) = ^'['X-i) = and w satisfies the second order linear ODE (13.4) on 
(xijXa). Therefore w = on (xi,X2), by uniqueness, which contradicts our 
construction of w. We conclude Ti < 0. D 

Motivation for the choice of trial function. The steady state equation reads 
U" + f (U) = 0, and differentiating shows that U' lies in the nullspace of the 
linearized operator L: 

LU' = -(U')"-f'(U)U' = 0. 

In other words, U' is an eigenfunction with eigenvalue 0, which almost proves 
instability (since instability would correspond to a negative eigenvalue). Of 
course, the eigenfunction U' does not satisfy the Dirichlet boundary condi- 
tions at the endpoints, and hence we must restrict to the subinterval (xi,X2), 
in the proof above, in order to obtain a valid trial function. 

Time maps and linearized stability 

Next we derive instability criteria that are almost necessary and sufficient. 
These conditions depend on the time map for a family of steady states. 

Parameterize the steady states by their slope at the left endpoint: given 
s y^ 0, write Ujx) for the steady state on M (if it exists) satisfying 

UJ0)=0, U;(0)=s, U^lx) =0 for some x>0. 

Define the time map to give the first point or "time" x at which the steady 
state hits the axis: 

T(s) = min{x > : Us(x) = 0}. 

If Us exists for some s 7^ then it exists for all nonzero s-values in a 
neighborhood, and the time map is smooth on that neighborhood [Schaaf 



Proposition 4.1.1]. The time map can be determined numerically by plotting 
solutions with different initial slopes, as the figures below show. In the first 
figure the time map is decreasing, whereas in the second it increases. 

Monotonicity of the time maps determines stability of the steady state: 
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C/,,+[/3=0 




Theorem 13.3 ( |Schaaf | Proposition 4.1.3]). The steady state Uj is lin- 
early unstable on the interval (0, T(s)) z/sT'(s) < 0, and is linearly stable if 
sT'(s) >0. 

Proof. We begin by differentiating the family of steady states with respect 
to the parameter s, and obtaining some properties of that function. Then 
we treat the "instabihty" and "stability" parts of the theorem separately. 

Write So 7^ for a specific value of s, in order to reduce notational con- 
fusion. Let X = T(so). Define a function 

9Us 



V 



9s 



S=So 



on (0,X), where we use that Us(x) is jointly smooth in (x, s). Then 

v" + f'(U)v = (13.5) 



as one sees by differentiating the steady state equation (13.2) with respect 
to s, and writing U for Usg. 
At the left endpoint we have 

v[0)=0, V(0) = 1, 

because Us(0) = 0,U;(0) = s for all s. 

We do not expect v to vanish at the right endpoint, but we can calculate 
its value there to be 

v(X) = sor(so), 

as follows. First, differentiating the equation = Us(T(s)) gives that 

^U3(T(s)) 

OS 

^^^(T(s)) + u;(T(s))r(s) 







9s 



v(T(s)) + u;(T(s))r(s). 
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Note the steady state Uj is symmetric about the midpoint of the interval 
(0,T(s)) (exercise; use that Us = at both endpoints and that the steady 
state equation is invariant under x h-> — x, so that steady states must be 
symmetric about any local maximum point). Thus 113(7(5)) = — Us(0) = 
— s, and evaluating the last displayed formula at s = Sq then gives that 
= v(X) — SoT'(so), as we wanted. 

Proof of instability. Assume SoT'(so) < 0. Thenv(X) < 0. Since v'(0) = 1 
we know v(x) is positive for small values of x, and so some Xz G (0, X) exists 
at which v(x2) = 0. Define a trial function 

J V on [0,X2), 
w = < 

I elsewhere. 

Then w is piecewise smooth, and is continuous since v = at X2. Note 
w(0) = 0. Therefore w e Hi(0,X), and w ^ 0. 



Hence Ti < by arguing as in the proof of Theorem |13.2[ except with 
xi =0. 

[Motivation for the choice of trial function. Differentiating the steady state 
equation U" + f (U) = with respect to s shows that 9U/9s is an eigenfunc- 
tion with eigenvalue zero: 

m m.._ m^^^ 

^ OS ^ OS OS 

In other words, 9U/9s lies in the nullspace of the linearized operator. It does 
not satisfy the Dirichlet boundary condition at the right endpoint, but we 
handled that issue in the proof above by restricting to the subinterval (0, X2), 
in order to obtain a valid trial function.] 

Proof of stability. Assume SoT'(so) > 0, so that v(X) > 0. Define ff = 
-v'(X)/v(X). Then 

v{0]=0, v'(X) + crv(X) = 0, 

which is a mixed Dirichlet-Robin boundary condition. We will show later 
that V is a first eigenfunction for L, under this mixed condition, with eigen- 



value is pi = (since Lv = by (13.5)) 



By adapting our Dirichlet-to- Robin monotonicity result (Theorem 10.1) 
one deduces that 

Ti > pi = 0, 
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which gives hnearized stabihty of the steady state U. 

To show V is a first eigenfunction for L, as used above, we first show v is 



positive on (0, X). Apply the steady state equation (13.2) to Us, and multiply 
by Ug and integrate to obtain the energy equation 



^-{u'f + nus 



1 



(13.6) 



0. Differentiating with 



where F is an antiderivative of f chosen with F(0) 
respect to s at s = Sq gives that 

UV + f (U)v = So. 



Hence if v vanishes at some Xo G [0, X) then UV — Sq ^ and so v'(xo) 7^ 
0. Thus at any two successive zeros of v, we know v' has opposite signs. 
Therefore U' has opposite signs too, because U'v' — so at the zeros. It is 



straightforward to show from ( 13.6 ) that U increases on [0, X/2] and decreases 



on [X/2, X] , and so after the zero of v at x = the next zero (if it exists) can 
only be > X/2, and the one after that must be > X. Since we know v(x) is 
positive for small x and that v(X) > 0, we conclude v has no zeros in (0, X) 
and hence is positive there. 

The first eigenfunction of L with mixed Dirichlet-Robin boundary con- 
dition is positive, and it is the unique positive eigenfunction (adapt the ar- 
gument in GilbargTrudinger , Theorem 8.38]). Since the eigenfunction v is 
positive, we conclude that it is the first Dirichlet-Robin eigenfunction, as 
desired. D 
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Chapter 14 

Case study: stability of steady 
states for thin fluid fllm PDEs 



Goal 

To linearize a particular nonlinear PDE around a steady state, and develop 
the spectral theory of the linearized operator. 



References LaugesenPughl , LaugesenPugh2 



Thin fluid film PDE 

The evolution of a thin layer of fluid (such as paint) on a fiat substrate (such 
as the ceiling) can be modeled using the thin fluid film PDE: 

Ht = -(f(h)hxxx)^ - (g(H)hx)^ 

where h,(x, t) > measures the thickness of the fluid, and the smooth, posi- 
tive coefficient functions f and g represent surface tension and gravitational 
effects (or substrate- fluid interactions). For simplicity we assume f = 1, so 
that the equation becomes 

Ht = -H^xxx - (g(H)hx)^. (14.1) 

We will treat the case of general g, but readers are welcome to focus on the 
special case g[y) =y^ for some p G M. 
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Solutions are known to exist for small time, given positive smooth initial 
data. But films can "rupture" in finite time, meaning h,(x, t) \^ as t ^ T, 
for some coefficient functions g (for example, for g = 1). 

Intuitively, the 4th order "surface tension" term in the PDE is stabilizing 
(since ht = — H^xxx is the usual 4th order diffusion equation) whereas the 
2nd order "gravity" term is destabilizing (since ht = — Hxx is the backwards 
heat equation). Thus the thin film PDE features a competition between 
stabilizing and destabilizing effects. This competition leads to nonconstant 
steady states, and interesting stability behavior. 

Periodic BCs and conservation of fluid. Fix X > and assume h, is 
X-periodic with respect to x. Then the total volume of fluid is conserved, 
since 



d 
dt 



pX rX 

h(x,t)dx = - {K^^^ + {g{h]K) J dx 
Jo 

x=X 

= -(h-xxx + g(H)Hx) 

x=0 



by periodicity. 

Nonconstant steady states. Every constant function is a steady state 



of (14.1). We discuss the stability of these steady states at the end of the 
chapter. 

To flnd nonconstant steady states, substitute h, — H(x) and solve: 

-Hxxxx-(g(H)Hx)x = (14.2) 

Hxxx + g(H)H, = a 
Hxx+G(H) = (3 + ax 
Hxx+G(H) = (3 

where G is an antiderivative of g; here a = is forced because the left side of 
the equation (Hxx + G(H)) is periodic. This last equation describes a nonlin- 
ear oscillator, and it is well known how to construct solutions (one multiplies 
by Hx and integrates). For example, when g = (27t/X)^ we have steady states 
H(x) = (const.) + cos(27tx/X). For the general case see |LaugesenPughl| . 
Assume from now on that H(x) is a nonconstant steady state with period X. 
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Linearized PDE 

We perturb a steady state by considering 

h= H + £c() 

where tlie perturbation 4)(x, t) is assumed to have mean value zer o J^ ^ {'X-, t) dx = 
0), so that fluid is conserved. Substituting h into the equation (14.1) gives 

0+£4)t = -(Hxxxx + e4^xxxx) - (g(H + £(t))(Hx + £4)x))^ 

= -Hxxxx - (g(H)H,)^ - £ [(t)xxx + g(H)4)x + g'(H)H,4)]^ + 0(£^). 

The leading terms, of order £°, equal zero by the steady state equation for 
H. We discard terms of order £^ and higher. The remaining terms, of order 
£\ give the linearized equation: 

*t = -[ct)xx + g(H)ct)]^^. (14.3) 

Unfortunately, the operator on the right side is not symmetric (meaning it 
does not equal its formal adjoint). To make it symmetric, we "integrate up" 
the equation, as follows. Write 

where i|) is X-periodic (since cj) has mean value zero). We may suppose ^\) 
has mean value zero at each time, by adding to \[> a suitable function of t. 
Substituting cj) = i|)x into ( 14. 3[ ) gives that 

A>tx = -[A>xxx + g(H)\|;J^^ 

^t = -[^xxx+ 9[^]^x]^ 

(noting the constant of integration must equal 0, by integrating both sides 
and using periodicity). Thus 

^\>, = -L^\> 

where L is the symmetric operator 

Lw = Wxxxx + (g(H)Wx)^. 
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Separation of variables gives (formally) solutions of the form 



^\) = ^qe ^''wjW) * = }! 



qe ^''w;(x), 



where the eigenvalues Tj and periodic eigenfunctions Wj satisfy 

LWj = TjWj. 

We conclude that the steady state H of the thin fluid film PDE is 

linearly unstable if Ti < 

because the perturbation 4) grows to infinity, whereas the steady state is 

linearly stable if Ti > 

because (f) remains bounded in that case. Remember these stability claims 
relate only to mean zero (volume preserving) perturbations. 

To make these claims more rigorous, we need to understand the eigenvalue 
problem for L. 

Spectrum of L 

We take: 

O = T = ]R/(XZ) = torus of length X, so that functions on D. are X- 
periodic 

"H = L^(T), inner product (u, v)l2 = Jq uv dx 

/C — H^(T) n {u : Jq u dx = 0}, with inner product 

pX 



(u,v)h2 = 



(u"v" + u'V + uv) dx 



Compact imbedding /C ^^ L^ by Rellich-Kondrachov 
Symmetric sesquilinear form 



a u,v 



f-X 





(uV - g(H)u'v' + Cuv) dx 

where C > is a sufficiently large constant to be chosen below. 
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Proof of ellipticity: The quantity a(u, u) has a term of the form — (u')^, 
whereas for ||u||^2 we need +(u')^. To get around this obstacle we "hide" 
the — (u')^ term inside the terms of the form (u")^ and u^. Specifically, 



( 





u'l^ dx 



rX 



u"u dx 



rX 



< 



(6(u")^ + (45)-^u^)dx 



(14.4) 



for any 5 > 0. Here we used "Cauchy-with-6, which is the observation that 
for any a, |3 G M, 



0< (V6a±(46)-^/^(3)^ 



1q2 



|a|3| < 5ot+ (46)"' (3 



Next, 
a(u,u) 

pX 

> 



((u- 



"\2 



|g(H)||Loo + l)(u')' + ^(u')' + Cu^) dx 



> 



([l-(||g(H)||Loc + l)5](u")' + l(u')'+[C-(||g(H)|koo + l)(45)->^)dx 



by (14.4) 



1,1 , 

> - u 



H2 



provided we choose 5 sufficiently small (depending on H) and then choose C 
sufficiently large. Thus ellipticity holds. 

The discrete spectral Theorem |4.1 now yields an ONB of eigenfunctions 
{wj} with eigenvalues yj such that 

a(wj,v) = yj(wj,v)i_2 Wv E JC. 

Writing yj = Tj + C we get 



rX pX 

(w;V'-g(H)w;v')dx = Tj 



WjV dx Vv G /C. 



These eigenfunctions satisfy Lwj = TjWj weakly, and hence also classically 
(by elliptic regularity, since H and g are smooth). 
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Zero eigenvalue due to translational symmetry. We will show that 
T = is always an eigenvalue, with eigenfunction u = H — H, where the 
constant H equals the mean value of the steady state H. Indeed, 

Lu = L(H-H) - H,,,,+ (g(H)H,)^ = 



by the steady state equation (14.2). 

This zero eigenvalue arises from a translational perturbation of the steady 
state, because choosing 

H = H(x + e) = H(x) + £H'(x) + 0(£^) 

gives rise to c() = H' and hence \[) = H — H. 



Sufficient condition for linearized instability of nonconstant steady 
state H 



Theorem 14.1 ( |LaugesenPugh2" Th. 3]). If g is strictly convex then r-\ < 0. 



For example, the theorem shows that nonconstant steady states are un- 
stable with respect to volume-preserving perturbations if gfy) = ij'' with 
either p > 1 or p < 0. 



Incidentally, the theorem is essentially the same as Theorem |13.1| for 
the reaction-diffusion PDE, simply writing g instead of f and noting that 
our periodic boundary conditions take care of the boundary terms in the 
integrations by parts. 

Proof. The Rayleigh principle for L says that 



Jo((w")'-g(H)(w')2)dx rx 



Ti = min <^ ^^^^ y — — : w G H^(T) \ {0} 



Jo ^^ d^ 
We choose 



wdx = 





w = H', 



which is not the zero function since H is nonconstant, and note w has mean 
value zero (as required), by periodicity of H. Then the numerator of the 
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Rayleigh quotient for w is 

'''((H'")'-g(H)(H")^)dx 



[■X 


rX 




1 pX 

3 



( - H"" - g ( H) H") H" dx by parts 

g'(H) (H')^H" dx by the steady state equation ( |14.2D 

g'(H)[(H')^'dx 
g"(H)(H'rdx 





3 



by parts 







<0 



by convexity of g and since H' ^ 0. Hence Ti < 0, by the Rayleigh principle. 

D 

Motivation for the choice of trial function. Our trial function w = H' corre- 
sponds to a perturbation cj) — H". This perturbation H = H + £(f) = H+ eH" 
tends to push the steady state towards the constant function. The oppo- 
site perturbation would tend to push the steady state towards a "droplet" 
solution that equals at some point. Thus our instability proof in Theo- 
rem 14.1 suggests (in the language of dynamical systems) that a heteroclinic 



connection might exist between the nonconstant steady state and the con- 
stant steady state, and similarly between the nonconstant steady state and 
a droplet steady state. 

Linear stability of nonconstant steady states. It is more difficult to 
prove stability results, because lower bounds on the first eigenvalue are more 
difficult to prove (generally) than upper bounds. 

See |LaugesenPugh2[ §3.2] for some results when gfy) — y'P, <j) < 3/4, 
based on time-map monotonicity ideas from the theory of reaction diffusion 



equations (see Chapter 13). 



Stability of constant steady states. 

Let H > be constant. Then H = H is a constant steady state. Its stability 
is easily determined, as follows. 
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Linearizing gives 

4>t = -^xxxx- g(H)(t)x 



by (14.3), where the right side is hnear and symmetric with constant coeffi- 
cients. 

We substitute the periodic Fourier mode (() = e^'^^exp(2mkx/X), where 
k G Z (and k ^ since our perturbations have mean value zero), obtaining 
the eigenvalue 

If g < (which means the second order term in the thin film PDE behaves 
like a forwards heat equation), then x > for each k, and so all constant 
steady states are linearly stable. 

If g > and (^) > g(H) then t > for each k, and so the constant 
steady states H is linearly stable. 

If g > and (y) < g(H) then the constant steady states H is linearly 
unstable with respect to the k = 1 mode (and possibly other modes too). In 
particular, this occurs if X is large enough. Hence we call the thin ffim PDE 
"long-wave unstable" if g > 0, since constant steady states are then unstable 
with respect to perturbations of sufficiently long wavelength. 
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Part II 



Continuous Spectrum 
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Looking ahead to continuous 
spectrum 



The discrete spectral theory in Part [T] of the course generated, in each apph- 
cation, 

• eigenfunctions {itj} with "discrete" spectrum Ai,A2,A3,... satisfying 
Litj = AjUj where L is a symmetric differential operator, together with 

• a spectral decomposition (or "resolution" ) of each f G L^ into a sum of 
eigenfunctions: f = ^j (f , ■U-j)Uj . 

These constructions depended heavily on symmetry of the differential opera- 
tor L (which ensured symmetry of the sesquilinear form a) and on compact- 
ness of the imbedding of the Hilbert space /C into l-i. 

For the remainder of the course we retain the symmetry assumption on 
the operator, but drop the compact imbedding assumption. The resulting 
"continuous" spectrum leads to a decomposition of f G L^ into an integral of 
"almost eigenfunctions" . 

We begin with examples, and later put the examples in context by de- 
veloping some general spectral theory for unbounded, selfadjoint differential 
operators. 
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Chapter 15 

Computable example: 
Laplacian (free Schrodinger) on 
all of space 



Goal 

To determine for the Laplacian on Euclidean space its continuous spectrum 
[0,oo), and the associated spectral decomposition of L^. 



Spectral decomposition 

The Laplacian L = —A on a bounded domain has discrete spectrum, as we 
saw in Chapters [2] and [5j When the domain expands to be all of space, 
though, the Laplacian has no eigenvalues at all. For example in 1 dimension, 
solutions of — u" = Au are linear combinations of e^^^"^", which oscillates 
if A > 0, or is constant if A = 0, or grows in one direction or the other if 
A G C \ [0,oo). In none of these situations does u belong to L^. (In all 
dimensions we can argue as follows: if — Au = Au and u G L^ then by taking 
Fourier transforms, 47t^|£,pu(£,) = u(£,) a.e., and so u = a.e. Thus no 
L^-eigenfunctions exist.) 

A fundamental difference between the whole space case and the case of 
bounded domains is that the imbedding H^ (M'^) "-^ L^(M'^) is not compact. 
For example, given any nonzero f G H^ (M), the functions f (■ — k) are bounded 
in L^(M), but have no L^-convergent subsequence as k ^ oo. Hence the 



discrete spectral theorem (Theorem 4.1) is inapplicable. 



Nevertheless, the Laplacian —A on M'^ has: 
1. generahzed eigenfunctions 



v„, X = e 



2ma)-x 



cu G 



(note that v^, is bounded, but it is not an eigenfunction since v^, ^ L^) 
which satisfy the eigenfunction equation — Avo, = Avcu with generahzed 
eigenvalue 



A = A(cu) =47r^| 



tu 



Generalized ei 




2. and a spectral decomposition 



(f,Va,)vtu dcu, Vf G L 



2fTn>d\ 



Proof of spectral decomposition. Since 



(f,Va,) 



f(x)e-""^"'-''dx = f(cu), 



the spectral decomposition simply says 



f(x) 



ftcuje^^'^-^dcu, 



which is the Fourier inversion formula. 
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Application of spectral decomposition 

One may solve evolution equations by separating variables: for example, the 
heat equation Ut = Au with initial condition h,(x) has solution 



u(x, t) 



h{w)e-^^'^^\a,M do). 



Note the analogy to the series solution by separation of variables, in the case 
of discrete spectrum. 

Aside. Typically, one evaluates the last integral (an inverse Fourier trans- 
form) and thus obtains a convolution of the initial data h, and the fundamen- 
tal solution of the heat equation (which is the inverse transform of e^'^''"'^). 

Continuous spectrum = [0, oo) 

The generalized eigenvalue A > is "almost" an eigenvalue, in two senses: 

• the eigenfunction equation [—A — A)u = does not have a solution in 
L^, but it does have a solution Vo, G L°°, 

• a Weyl sequence exists for —A and A, meaning there exist functions 
Wn such that 

(Wl) ||(— A — A)wn||L2 ^ as n ^ CO, 

(W2) ||Wn||L2=1, 

(W3) Wn —^ weakly in L^ as n ^ oo. 

We prove existence of a Weyl sequence in the Proposition below. Later we 
will define the continuous spectrum to consist of those A-values for which 
a Weyl sequence exists. Thus the continuous spectrum of —A is precisely 
the nonnegative real axis. Recall it is those values of A that entered into our 
spectral decomposition earlier in the chapter. 

Remark. Existence of a Weyl sequence ensures that (—A — A) does not have 
a bounded inverse from L^ -^ L^, for if we write fn = (—A — A)Wn then 

||(-A-A)-^fn||L2 l|Wn||L2 
^ — } CX) 

||fn||L2 IK-A-A)Wn||L2 

as TL ^ OO, by (Wl) and (W2). In this way, existence of a Weyl sequence is 
similar to existence of an eigenfunction, which also prevents invertibility of 
(-A-A). 
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c 



spec(-A) 



Proposition 15.1 (Weyl sequences for negative Laplacian). A Weyl se- 
quence exists for —A and A G C z/ and only if7\& [0, oo). 



Fix A G [0, cxd) and choose cu G M'^ witli 47t|a)p = A. Take a 
'^] sucii that K = 1 on the unit ball B(1 ) and k = 
'^ \ B(2). Define a cut-off version of the generahzed eigenfunction, by 



Proof. 

cut-off function k G C^ 

on ™-' ' "'^" " " 



w^ 



CnK 



v„ 



n 



where the normalizing constant is 



n 



d/21 



K||l2 



Weyl approx. eigenfn. W4(x) 




First we prove (Wl). We have 

(A + A)Wn 
= Cn(AVa, + AVa,)K(-) +2^VV^(X) • (Vk)(-) + ^V^ (x) (Ak) ' "^ ' 



n 



n^ 



n 



n n 

The first term vanishes because Avw = — 47t|cupVa, pointwise. In the third 
term, note that Vo, is a bounded function, and that a change of variable shows 



n^ n TL^ ||k||l2 



0. 
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The second term similarly vanishes in the limit, as n -^ oo. Hence (A + 
A]wn ^ in L^, which is (Wl). 

For (W2) we simply observe that |Va)(x)| = 1 pointwise, so that ||wn.||L2 = 
1 by a change of variable, using the definition of c^- 

To prove (W3), take f G L^ and let R > 0. We decompose f into "near" 
and "far" components, as f = g + H where g = flB(R) and h, = f1]Rd\B(R). 
Then 

(f,Wn)L2 = (g,Wn)L2 + {h,Wn)l2. 

We have 

|(g,Wn)L2| < Cn||K||L-||g||L' -^ 

as TL ^^ oo, since Cn -^ 0. Also, by Cauchy-Schwarz and (W2) we see 

hmsup |(h,,Wn)L2| < ||H||l2. 

n— >oo 

This last quantity can be made arbitrarily small by letting R ^ oo, and so 
limn-_>oo(f) Wn)L2 = 0. That is, Wn ^ weakly. 

"=>" Assume A G C \ [0, oo), and let 

5 = dist (A, [0,oo)) 

so that 6 > 0. 

Suppose (Wl) holds, and write gn = (— A — A)wn. Then 

g;(£,) = (47t2|£,|2-A)T^(£,) 
-ll 



and hence 



\^n{^)\<&-'\g^{^)\ 



|Wn||L2 = ||w^||l2 < s lg^||L2 



= 5 lgn||L2 

by (Wl). Thus (W2) does not hold. 

{Aside. The calculations above show, in fact, that (— A — A)^^ is bounded 
from L^ -^ L^ with norm bound 8^\ when A ^ [0, oo).) D 
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Chapter 16 

Computable example: 
Schrodinger with a bounded 
potential well 



Goal 

To show that the Schrodinger operator 



^ — 2 sech X 



in 1 dimension has a single negative eigenvalue (discrete spectrum) as well 
as nonnegative continuous spectrum [0, oo). The spectral decomposition will 
show the potential is refiectionless. 



-2 sech X potential 




C 



-1 



spe<;(L) 



Reference |Keenerj Section 7.5 
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Discrete spectrum = {—1} 

We claim —1 is an eigenvalue of L with eigenfunction sechx. This fact can 
be checked directly, but we will proceed more systematically by factoring the 
Schrodinger operator with the help of the first order operators 

T + d 1 

L = h tanhx, 

dx 

L^ = - — h tanh x. 
dx 



We compute 



L+L" - 1 = (--— + tanhx) (-— + tanhx) - 1 
^ dx ^ Mx ' 

d^ 

(tanh x)' + tanh^ x — 1 



dx2 

_d^ 

dx2 
L 



2 sech^ X 



since (tanh)' — sech and 1 — tanh = sech . Thus 

L = L+L--1. (16.1) 

It follows that functions in the kernel of L^ are eigenfunctions of L with 
eigenvalue A = — 1 . To find the kernel we solve: 

v' + (tanhx)v = 
(coshx)v' + (sinhx)v = 

(coshx)v = const. 
V = c sech X 

Clearly sechx G l-^(M), since sech decays exponentially. Thus —1 lies in the 
discrete spectrum of L, with eigenfunction sechx. 

Are there any other eigenvalues? No! Argue as follows. By composing 
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Bound state, with energy -1 




L+ and L in the reverse order we find 



L L^ 



( h tanh x) ( h tanh x) — 1 

Mx ^ ^ dx ' 

Al 

+ (tanhx)' + tanh^ x — 1 



dx2 
d^ 



dx^' 



:i6.2) 



Prom (16.1) and (16.2) we deduce 



dx^ 



L L. 



Thus if Lv = Av then -^ 



L v) = L Lv = A(L v). By solving for L v in 

terms of e^^^"^", and then integrating to obtain v, we conclude after some 
thought (omitted) that the only way for v to belong to L^(]R) is to have 



L V = and hence v — c sech x, so that A 



-1. 



Continuous spectrum D [0, oo) 

Let A G [0,oo). Generalized eigenfunctions with Lv = Av certainly exist: 
choose cu G M with A = 47r^a)^ and define 



vfxl = L+fe' 



(tanhx — 2niw)e^ 
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which is bounded but not square integrable. We compute 



Lv = (L+L- - IjL+le^"'^") by ( |16.lD 

= L+(L-L+-l)(e^"^^'') 
d^ 



-L^^le^"^-") by (|163 

= L+ (471^0)^6^™^") 

= Av, 

which verifies that v(x) is a generahzed eigenfunction. 

We can further prove existence of a Weyl sequence for L and A by adapting 



Lemma 15.1 "<^", using the same Weyl functions Wn(x) as for the free 
Schrodinger operator —A. The only new step in the proof, for proving ||(L — 
A)wtt^||]_2 — > in (Wl), is to observe that 

|2sech^xWn(x)|=2Cn|K(-)e^"^'"''|sech2x 

n 

< 2Cn||K||i_cx> sech X 

in L^(]R) as n ^ oo, because Cn -^ 0. {Note. This part of the proof works 
not only for the sech potential, but for any potential belonging to L^.) 

We have shown that the continuous spectrum contains [0, oo). We will 
prove the reverse inclusion at the end of the chapter. 

Generalized eigenfunctions as traveling waves. The eigenfunction 
("bound state") v(x) = sechx with eigenvalue ("energy") —1 produces a 
standing wavefunction 

u = e^^ sech x 

satisfying the time-dependent Schrodinger equation 

iut = Lu. 
The generalized eigenfunction 

v(x) = (tanhx-2mcu)e^"^^'' (16.3) 

with generalized eigenvalue A = 47r^uj^ similarly produces a standing wave 
u = e-^'*"'"'''(tanhx - 2ma))e^"^"'\ 
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Gen. eigenfn, ^= — : Re v(x) 




More usefully, we rewrite this formula as a traveling plane wave multiplied 
by an x-dependent amplitude: 



u = (tanhx - 2ma;)e^™'"t''-^"^''. 



(16.4) 



The amplitude factor serves to quantify the effect of the potential on the 
traveling wave: in the absence of a potential, the amplitude would be identi- 
cally 1 , since the plane wave e^Tticutx-iTtajt) gQ^ygg ^^q fj-gg Schrodinger equation 
iitt = — Au. 

Refiectionless nature of the potential, and a nod to scattering the- 
ory. One calls the potential — 2sech x "refiectionless" because the right- 



moving wave in (16.4) passes through the potential with none of its energy 



reflected into a left-moving wave. In other words, the generalized eigenfunc- 
tion (16.3) has the form ce^"^'^" both as x ^ — oo and as x ^ oo (with 



different constants, it turns out, although the constants are equal in magni- 
tude). 

This refiectionless property is unusual. A typical Schrodinger potential 
would produce generalized eigenfunctions equalling approximately 



j,^g2ma,x ^ c^Q-^nia^x 



as X 



-CXD 



and 



CtB 



Inxwx 



as X 



oo 



(or similarly with the roles of ±oo interchanged). Here |ci| is the amplitude 
of the incident right-moving wave, |cr| is the amplitude of the left-moving 
wave reflected by the potential, and |ct| is the amplitude of the right-moving 
wave transmitted through the potential. Conservation of L^-energy demands 
that 

|ci|' = |cr|2 + |ct|'. 
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For a gentle introduction to this "scattering theory" see [Keener! Section 7.5]. 
Then one can proceed to the book-length treatment in |ReedSinion3j . 

Spectral decomposition of L^ 

Analogous to an orthonormal expansion in terms of eigenfunctions, we have: 



Theorem 16.1. 



1 



r 



2 



dcu 



f = :r(f,sech) sech+ (f, L+Va,)L+v,, ^^^, Vf G L 



1 +4n^uo^' 



where L+V|^(x) = (tanhx — Imcuje^"^^" is the generalized eigenfunction at 
frequency tu. 

The discrete part of the decomposition has the same form as the contin- 
uous part, in fact, because sech = — L+(sinh). 

Proof. We will sketch the main idea of the proof, and leave it to the reader 
to make the argument rigorous. 

By analogy with an orthonormal expansion in the discrete case, we assume 
that f G L^(]R) has a decomposition in terms of the eigenfunction sechx and 
the generalized eigenfunctions L+v^, in the form 

f = c(f, sech) sech+ raf(a))(f, L^Va,)L+Vcu dcu, 

Jr 

where the coefficient c and multiplier raf(cu) are to be determined. 

Taking the inner product with sech x implies that c = |, since ||sech||^2f]gi = 
2 and (L+Va,,sech) = (va,,L^sech) = 0. 

Next we annihilate the sech term by applying L^ to both sides: 

L-f = L-( mf(a))(f,L+v^)L+v^dcu 
^ Jr 

Note that by integration by parts, 

(f,L+v^) = (L-f,Va,) = rL^)(cu). 
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Hence 



L f = L 

r 



raf(cu)(L f)(uj)L v^, dcu 
raf(uj)(L"f)(uj)L^L+Va, dcu 



by (16.2). Thus the multipher should be raf(uj) = 1/(1 +4n^(v^], in order 



for Fourier inversion to hold. This argument shows the necessity of the 
formula in the theorem, and one can show sufficiency by suitably reversing 
the steps. D 

The theorem implies a Plancherel type identity. 

Corollary 16.2. 

dw 



f||f, = ^Kf.Bech)|^ + ||(f,L-v„)P^^^^, 



Vf G L^( 



Proof. Take the inner product of f with the formula in Theorem 16.1 D 



Continuous spectrum = [0, oo) 

Earlier we showed that the continuous spectrum contains [0, oo). For the 
reverse containment, suppose A ^ [0, oo) and A ^ — 1 . Then L— A is invertible 
on L^, with 



:L-A)-'f 



1 1 



A+12 



(f , sech) sech + 



f- (f,L+v^) ^ dcv 

4n^w^ - A ^"^ 1 + 4n^cv^ 



as one sees by applying L — A to both sides and recalling Theorem 16.1 To 
check the boundedness of this inverse, note that 



[L-A)-'f 



1-CII2 



'^' IA+1P2 

1 1 
< 



|(f,sech)|^ + 



|(f,L+v^)P dcu 



|A+1|22 
< (const.)||f ||]_2, 



|(f,sech)p + 



\4n^w^ - AM + 4n^cv^ 



dist (A, [0,oo))^ 



dcu 



1 + 4n^w^ 
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where we used Corollary |16.2 



The boundedness of (L — A)^^ implies that the Weyl conditions (Wl) and 
(W2) cannot both hold. Thus no Weyl sequence can exist for A, so that A 
does not belong to the continuous spectrum. 

Next suppose A = — 1 . If a Weyl sequence w^ exists, then 

(wn, sech)i_2 -^ as n ^ oo, 

by the weak convergence in (W3). Hence if we project away from the A = — 1 
eigenspace by defining 

Un = Wn - ^('^^n, sech)L2 sech and z^ = Un/llynllL^) 

then we find HUtiIIl^ -^ 1 and ||zn.||L2 = 1, with (zn,sech)]_2 = 0. Also 

(L+ 1)Zn = (L+ l)lJn/||yn||L2 = (L+ 1)Wn/||yn||L2 ^ 

in L^. Thus Zn satisfies (Wl) and (W2) and lies in the orthogonal com- 
plement of the eigenspace spanned by sech. A contradiction now follows 
from the boundedness of (L+ 1)^^ on that orthogonal complement (with the 
boundedness being proved by the same argument as above for A 7^ —1 ). This 
contradiction shows that no such Weyl sequence w^ can exist, and so —1 
does not belong to the continuous spectrum. 

Note. The parallels with our derivation of the continuous spectrum for 



the Laplacian in Chapter 15 are instructive. 
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Chapter 17 

Selfadjoint, unbounded linear 
operators 



Goal 

To develop the theory of unbounded linear operators on a Hilbert space, and 
to define selfadjointness for such operators. 



References |GustafsonSigal| Sections 1.5, 2.4 



^HislopSigalj Chapters 4, 5 



Motivation 

Now we should develop some general theory, to provide context for the ex- 



amples computed in Chapters 15 and 16 



We begin with a basic principle of calculus: 

integration makes functions better, while differentiation makes 
them worse. 

More precisely, integral operators are bounded (generally speaking), while 
differential operators are unbounded. For example, e^"^^^" has norm 1 in 
L^[0, 1] while its derivative i-e^"""^" — Zmne^™"^" has norm that grows with 
n. The unboundedness of such operators prevents us from applying the 
spectral theory of bounded operators on a Hilbert space. 

Further, differential operators are usually defined only on a (dense) sub- 
space of our natural function spaces. In particular, we saw in our study 
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of discrete spectra that the Laplacian is most naturally studied using the 
Sobolev space H\ even though the Laplacian involves two derivatives and 
H^ -functions are guaranteed only to possess a single derivative. 

To meet these challenges, we will develop the theory of densely defined, 
unbounded linear operators, along with the notion of adjoints and selfad- 
jointness for such operators. 

Domains and inverses of (unbounded) operators 

Take a complex Hilbert space "H with inner product (■,■). Suppose A is a 
linear operator (not necessarily bounded) from a subspace D(A) C "H into 
H: 

A:D(A) ^n. 

Call D(A) the domain of A. 

An operator B with domain D(B) is called the inverse of A if 

• D(B) =Ran(A),D(A) =Ran(B), and 

• BA = idRan(B)) AB = idRan(A)- 

Write A^^ for this inverse, if it exists. Obviously A^^ is unique, if it exists, 
because in that case A is bijective. 

Further say A is invertible if A^^ exists and is bounded on "H (mean- 
ing that A^^ exists, Ran(A) = "H, and A^^ : "H ^ "H is a bounded linear 
operator). 

Example. Consider the operator A = — A+ 1 with domain H^(]R'^) C L^(M'^). 
Invertibility is proved using the Fourier transform: let D(B) = L^(M'^), and 
define a bounded operator B : L^ ^ L^ by 

One can check that Ran(B) = H^(M'i) = D(A). Notice BA = idH2,AB = 
idL2. The second identity implies that Ran(A) — \J- — D(B). 

Adjoint of an (unbounded) operator 
Call A symmetric if 

(Af,g) = (f,Ag), Vf,gGD(A). (17.1) 
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Symmetry is a simpler concept than selfadjointness, which requires the 
operator and its adjoint to have the same domain, as we now explain. 
First we define a subspace 

D(A*) = {f G "H : the linear functional g ^-^ (f) Ag) is bounded on D(A)}. 

Assume from now on that A is densely defined, meaning D(A) is dense in 
"H. Then for each f G D(A*), the bounded linear functional g i-^ (f, Ag) is 
defined on a dense subspace of Ti and hence extends uniquely to a bounded 
linear functional on all of "H. By the Riesz Representation Theorem, that 
linear functional can be represented as the inner product of g against a unique 
element of T-L, which we call A*f . Hence 

(f,Ag) = (A*f,g), Vf G D(A*), g G D(A). (17.2) 

Clearly this operator A* : D(A*) ^ "H is hnear. We call it the adjoint of A. 

Lemma 17.1. If A is a densely defined linear operator and A G C, then 
(A -A)* = A*-A. 

We leave the (easy) proof to the reader. Implicit in the proof is that 
domains are unchanged by subtracting a constant: D(A — A) = D(A) and 
D(((A-A)*) =D(A*). 

The kernel of the adjoint complements the range of the original operator, 
as follows. 



Proposition 17.2. If A is a densely defined linear operator then Ran(A) © 
ker(A*) =n. 

Proof. Clearly ker(A*) C Ran(A)-^, because if f G ker(A*) then A*f — and 
so for all g G D(A) we have 

(f,Ag) = (A*f,g)=0. 

To prove the reverse inclusion, RanfA)-*- C ker(A*), suppose h, G RanfA)-*-. 
For all g G D(A) we have (h, Ag) = 0. In particular, h, G D(A*). Hence 

(A*h,g) = (h,Ag)=0 Vg G D(A), 

and so from density of D(A) we conclude A*h, = 0. That is, H G ker(A*). 

We have shown RanfA)-*" = ker(A*), and so (since the orthogonal com- 
plement is unaffected by taking the closure) Ran(A) = ker(A*). The propo- 
sition follows immediately. D 
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We will need later that the graph of the adjoint, {(f, A*f) : f G D(A*)}, 
is closed in Ti x Ti. 

Theorem 17.3. If A is a densely defined linear operator then A* is a closed 
operator. 

Proof. Suppose fn G D(A*) with fn -^ f, A*fn -^ g, for some f , g G "H. To 
prove the graph of A* is closed, we must show f G D(A*) with A*f = g. 
For each h G D(A) we have 

(f, Ah) = lim(fn, Ah) = lim(A*fn,h) = (g,h). 

n n 

Thus the map h \-^ (f, Ah) is bounded for h G D(A). Hence f G D(A*), and 
using the last calculation we see 

(A*f,h) = (f,Ah) = (g,h) 

for all h G D(A). Density of the domain implies A*f = g, as we wanted. D 

Selfadjointness 

Call A selfadjoint if A* = A, meaning D(A*) = D(A) and A* = A on their 
common domain. 

Selfadjoint operators have closed graphs, due to closedness of the adjoint 
in Theorem 117.31 Thus: 

Proposition 17.4. // a densely defined linear operator A is selfadjoint then 
it is closed. 

The relation between selftadjointness and symmetry is clear: 

Proposition 17.5. The densely defined linear operator A is selfadjoint if 
and only if it is symmetric and D(A) = D(A*). 



Proof. "^=^" If A* = A then the adjoint relation (17.2) reduces immediately 



to the symmetry relation (17.1 ). 



The symmetry relation (17.1) together with the adjoint relation 



(17.2) implies that (Af , g) = (A*f, g) for all f , g G D(A) = D(A*). Since 



D(A) is dense in H, we conclude Af = A*f. D 

For bounded operators, selfadjointness and symmetry are equivalent. 
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Lemma 17.6. // a linear operator A is bounded on T-i, then it is selfadjoint 
if and only if it is symmetric. 

Proof. Boundedness of A ensures that D(A*) = "H = D(A), and so the 



adjoint relation (17.2) holds for all f, g G "H. Thus A* = A is equivalent to 



symmetry. D 

Example: selfadjointness for Schrodinger operators 

Let L — —A + V be a Schrodinger operator with potential V(x) that is 
bounded and real- valued. Choose the domain to be D(L) = H^(]R'^) in the 
Hilbert space L^(]R'^). This Schrodinger operator is selfadjoint. 

Proof. Density of D(L) follows from density in L^ of the smooth functions 
with compact support. 

Our main task is to determine the domain of L*. Fix f , g G H^(]R'^). 
From the integration by parts formula (f, Ag)i_2 = (Af, g)i_2 (which one may 
alternatively prove with the help of the Fourier transform) , one deduces that 

|(f,Ag)L2| = |(Af,g)L2|< ||f||H2||g||L2. 

Also |(f, Vg)i_2| < ||f ||]_2||V||i_cx,||g||i_2. Hence the linear functional g i— > (f,Lg)i_2 
is bounded on g G D(L). Therefore f G D(L*), which tells us H^(]R'^) C 
D(L*). 

To prove the reverse inclusion, fix f G D(L*). Then 

Kf,Lg)L2| < (const.)||g||L2, Vg G D(L) = H^(M'i). 

Since the potential V is bounded, the last formula still holds if we replace V 
with 1 , so that 

Kf, (-A+ l)g)L2| < (const.)||g||L2, Vg G H^fM'^). 

Taking Fourier transforms gives 

|(f,(1 +47t^|£,hg)L2| < (const.)||g||L2, Vg G H'IR"]. 

In particular, we may suppose g = h G C^(]R'^), since every such g gives 
g G H^tM'i). Hence 

|((1 +47r2|£,|^)f,h)L2| < (const.) II h||L2, Vhe C^(M^). 
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Taking the supremum of the left side over all h. with L^-norm equal to 1 
shows that 

11(1 +47t^|£,|^)^|L2 < (const.) 

Hence [^+\E,\)^?e L2(M'1), which means f G H2(M'1). Thus D(L*) C H^[R'^). 
Now that we know the domains of L and L* agree, we have only to check 
symmetry, and that is straightforward. When f , g G H^(]R'^) we have 

(Lf,g)=-(Af,g)L2 + (Vf,g)L2 

= -(f,Ag)L2 + (f,Vg)L2 

where we integrated by parts and used that V(x) is real-valued. D 
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Chapter 18 

Spectra: discrete and 
continuous 



Goal 

To develop the spectral theory of selfadjoint unbounded linear operators. 



References [GustafsonSigal Sections 2.4, 5.1 



jHislopSigall Chapters 1, 5, 7 
|Rudinj Chapter 13 

Resolvent set, and spectrum 

Let A be a densely defined linear operator on a complex Hilbert space "H, 
as in the preceding chapter. The operator A — A has domain D(A), for each 
constant A G C Define the resolvent set 

res(A) = {A G C : A — A is invertible (has a bounded inverse defined on T-L)}. 

For A in the resolvent set, we call the inverse (A — A)^^ the resolvent op- 
erator. 

The spectrum is defined as the complement of the resolvent set: 

spec(A) =C\res(A). 

For example, if A is an eigenvalue of A then A G spec (A), because if Af = Af 
for some f 7^ 0, then (A — A)f = and so A — A is not injective, and hence 
is not invertible. 
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Proposition 18.1 ( [HislopSigal Theorem 1.2]). The resolvent set is open, 



and hence the spectrum is closed. 

We omit the proof. 

The next result generahzes the fact that Hermitian matrices have only 
real eigenvalues. 

Theorem 18.2. If A is selfadjoint then its spectrum is real: spec(A) C M. 

Proof. We prove the contrapositive. Suppose A G C has nonzero imaginary 
part, ImA y^ 0. We will show A G res(A). 

The first step is to show A — A is injective. For all f G D(A), 



and so 



(A-A)f|^ = \\Af\^- 


-2(ReA)(f,Af) + |Anf 1^ 


(A-A)f |2 > |Af |^- 
=(|Af|- 
> llmAl^l 


-2|ReA| f 1 Af | + |A|2 f ' 
-|ReA||f|)^ + |ImA|^ f |^ 



(18.1) 

The last inequality implies that A — A is injective, using here that | Ini A| > 0. 

That is, ker( A - A) = {0}. 

Selfadjointness (A* = A) now gives kerf A*— A) = 0, and so Ran(A — A) = 



Ti by Proposition |17.2[ That is, A — A has dense range. 

Next we show Ran(A — A) = "H. Let g G "H. By density of the range, we 
may take a sequence fn G D(A) such that (A — Ajf^x -^ g. The sequence f^ 



is Cauchy, in view of (18.1). Hence the sequence (fn, (A — A)fn) is Cauchy 
in Ti X Ti, and so converges to (f, g) for some f G "H. Note each ordered 
pair (fn, (A — A)fn) lies in the graph of A — A, and this graph is closed 



by Proposition 17. 4| (relying here on selfadjointness again). Therefore (f, g) 



belongs to the graph of A — A, and so g G Ran(A — A). Thus A — A has full 
range. 

To summarize: we have shown A — A is injective and surjective, and so it 
has an inverse operator 

(A-A)-i :H^D(A) cn. 

This inverse is bounded with 

||(A-A)-ig||<|lmArlg||, WgeU, 



by taking f = (A — A) g in estimate (18.1). The proof is thus complete. D 
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Characterizing the spectrum 

We will characterize the spectrum in terms of approximate eigenfunctions. 
Given a number A G C and a sequence Wn G D(A), consider three conditions: 

(Wl) ||(A — A)wn||^ — > as n ^ CO, 

(W2) ||Wn||w = 1, 

(W3) Wn -^ weakly in "H as n ^ cxd. 



(We considered these conditions in Chapter 15 for the special case of the 
Laplacian). 

Condition (Wl) says Wn is an "approximate eigenfunction" , and condi- 
tion (W2) simply normalizes the sequence. These conditions characterize the 
spectrum, for a selfadjoint operator. 

Theorem 18.3. If A is selfadjoint then 

spec(A) = {A G C : (Wl) and (W2) hold for some sequence Wn G D(A)}. 

Proof, "d" Assume (Wl) and (W2) hold for A, and that A — A has an inverse 
defined on Ti. Then for f^ = (A — Ajw^^ we find 

||(A-A)"^fn||^ \\yvn\\n 

' oo 



||fn|k IKA-A)Wn||^ 

as n ^ cxD, by (Wl) and (W2). Thus the inverse operator is not bounded, 
and so A G spec(A). 



"c" Assume A G spec(A), so that A is real by Theorem 18.2, If A is an 
eigenvalue, say with normalized eigenvector f , then we simply choose w^ = f 
for each n, and (Wl) and (W2) hold trivially. 

Suppose A is not an eigenvalue. Then A — A is injective, hence so is 
(A — A)*, which equals A — A by selfadjointness of A and reality of A. Thus 



17.2 



ker ((A — A)*) = {0}, and so Ran(A — A) is dense in "H by Proposition 

Injectivity ensures that (A— A)^^ exists on Ran(A— A). If it is unbounded 
there, then we may choose a sequence fn G Ran(A— A) with || (A— A)^^fn||^ = 
1 and \\fn\\n -^ 0. Letting w^ = (A— A)^^fn gives (Wl) and (W2) as desired. 
Suppose on the other hand that (A — A)^^ is bounded on Ran(A — A). Then 



the argument in the proof of Theorem 18.2 shows that Ran(A — A) = "H, 
which means A belongs to the resolvent set, and not the spectrum. Thus this 
case cannot occur. D 
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Discrete and continuous spectra 

Define the discrete spectrum 

= {A G spec (A) : A is an isolated eigenvalue of A having finite multiplicity}, 

where "isolated" means that some neighborhood of A in the complex plane 
intersects spec (A) only at A. By "multiplicity" we mean the geometric mul- 
tiplicity (dimension of the eigenspace); if A is not selfadjoint then we should 
use instead the algebraic multiplicity |HislopSigal| . 
Next define the continuous spectrum 

SPeCcontl'^) 

= {A G C : (Wl), (W2) and (W3) hold for some sequence Wn G D(A)}. 



The continuous spectrum lies within the spectrum, by Theorem 18.3 The 
characterization in that theorem required only (Wl) and (W2), whereas the 
continuous spectrum imposes in addition the "weak convergence" condition 
(W3). 

A Weyl sequence for A and A is a sequence w^ G D(A) such that (Wl), 
(W2) and (W3) hold. Thus the preceding definition says the continuous 
spectrum consists of A- values for which Weyl sequences exist. 

The continuous spectrum can contain eigenvalues that are not isolated 
( "imbedded eigenvalues" ) or which have infinite multiplicity. 

A famous theorem of Weyl says that for selfadjoint operators, the entire 
spectrum is covered by the discrete and continuous spectra. 

Theorem 18.4. If A is selfadjoint then 

spec (A) = spec^^^jA) U spec,„^jA). 

(Further, the discrete and continuous spectra are disjoint.) 

We omit the proof. See [HislopSigal , Theorem 7.2]. 



Applications to Schrodinger operators 

The continuous spectrum of the Laplacian —A equals [0,oo), and the spec- 
trum contains no eigenvalues, as we saw in Chapter [15J 
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The hydrogen atom too has continuous spectrum [0, oo), with its Schrodinger 
operator L = — A — 2/|x| on M? having domain H^(M'^) C L^(]R'^); see [Taylor 



Section 8.7]. The discrete spectrum {— l/n^ : n > 1} of the hydrogen atom 
was stated in Chapter [3] 

As the hydrogen atom example suggests, potentials vanishing at infinity 
generate continuous spectrum that includes all nonnegative numbers: 

Theorem 18.5. Assume V[x] is real-valued, continuous, and vanishes at 
infinity (V[x] -^ as \x\ — > ooj. 

Then the Schrodinger operator —A+V is self adjoint (with domain H^(]R'^) C 
L^(M'^)j and has continuous spectrum = [0,oo). 



For a proof see [HislopSigal Corollary 14.10], where a stronger theorem 



is proved that covers also the Coulomb potential — 2/|x| for the hydrogen 
atom. Note the Coulomb potential vanishes at infinity but is discontinuous 
at the origin, where it blows up. The stronger version of the theorem re- 
quires (instead of continuity and vanishing at infinity) that for each e > 0, 
the potential V(x) be decomposable as V = V2 + Voo where V2 G L^ and 
||Voo||l°° < £• This decomposition can easily be verified for the Coulomb 
potential, by "cutting off" the potential near infinity. 

Theorem 18. 5| implies that any isolated eigenvalues of L must lie on the 



negative real axis (possibly accumulating at 0). For example, the — 2sech 
potential in Chapter [16] generates a negative eigenvalue at —1. 



Connection to generalized eigenvalues and eigenfunctions 

Just as the discrete spectrum is characterized by eigenfunctions in L^, so the 
full spectrum is characterized by existence of a generalized eigenfunction that 
grows at most polynomially at infinity. 

Theorem 18.6. Assume V(x) is real-valued and bounded on M.'^. Then the 
Schrodinger operator —A + V has spectrum 

spec(-A + V) = 

closure o/{A G C : (—A + V)u = Au for some polynomially bounded u}. 



We omit the proof; see GustafsonSigal Theorem 5.22]. 

Ill 



Further reading 

A wealth of information on spectral theory, especially for Schrodinger oper- 
ators, can be found in the books |GustafsonSigal[ |HislopSigai] IReedSimon2[ 
IReedSimon4] . 
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Chapter 19 

Discrete spectrum revisited 



Goal 



To fit the discrete spectral Theorem 4.1 (from Part ll] of the course) into 



the spectral theory of selfadjoint operators and, in particular, to prove the 
absence of continuous spectrum in that situation. 

Discrete spectral theorem 



The discrete spectral Theorem |4.1| concerns a symmetric, elliptic, bounded 
sesquilinear form a(u, v) on an infinite dimensional Hilbert space /C, where 
/C imbeds compactly and densely into the Hilbert space %. The theorem 
guarantees existence of an ONB for T-i consisting of eigenvectors of a: 

a(uj,v) =yj(uj,v)^ Vv G /C, 

where the eigenvalues satisfy 

< yi < 72 < 73 < ■ ■ ■ ^ oo- 

We want to interpret these eigenvalues as the discrete spectrum of some 
selfadjoint, densely defined linear operator on "H. By doing so, we will link 
the discrete spectral theory in Part |T] of the course with the spectral theory 
of unbounded operators in Part [TTj 

Our tasks are to identify the operator A and its domain, to prove A is 
symmetric, to determine the domain of the adjoint, to conclude selfadjoint- 
ness, and finally to show that the spectrum of A consists precisely of the 
eigenvalues yj. 

113 



Operator A and its domain 



In the proof of Theorem 4.1 we found a bounded, selfadjoint hnear operator 



B : H ^ /C C H with eigenvalues 1/yj and eigenvectors Uj: 

1 



Buj = — Uj. 



We showed B is injective (meaning its eigenvalues are nonzero). Notice B 
has dense range because its eigenvectors xtj span l-L. 

{Aside. This operator B relates to the sesquilinear form a by satisfying 
a(Bf, v) = (f, v)^ for all v G /C. We will not need that formula below.) 

Define 

A = B"^ Ran(B) -^ H. 

Then A is a linear operator, and its domain 

D(A) =Ran(B) 
is dense in %. 

Symmetry of A 

Let u,v G D(A). Then 

(Au,v)^ = (Au, BAv)^ since BA = Id, 

= (BAu, Av)^ since B is selfadjoint, 

= (u, Av)^ since BA = Id. 

Domain of the adjoint 

First we show D(A) C D(A*). Let u G D(A). For all v G D(A) we have 

I (u, Av) ^ I = I ( Au, v)u\ by symmetry 

< ||Au||^||v||^. 

Hence the functional v h-> (u, Av)^ is bounded on D(A) with respect to the 
?{-norm, so that u belongs to the domain of the adjoint A*. 
Next we show D(A*) C D(A). Let u G D(A*) C ?^. We have 

|(u,Av)>^| < (const.) ||v||^ Vv G D(A) = Ran(B). 
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Writing v = Bg gives 

l(u,g)Hl < (const. )||Bg||^ Vg G "H. 

One can express u in terms of the ONB as u = ^- djUj. Fix J > 1 and 

choose g = ^Li yf djUj G "H, so that Bg = Y.l=]'Yi^]'^]- We deduce from 
the last inequahty that 

^yf|d^|^<(const.)(^yf|dJp)^/^ 

and so 

J 

^yfldjP < (const. )^ 

j=i 

Letting J ^ oo imphes that 

^yfldjl^ < (const. )^ 
j 

and so the sequence {yjdj} belongs to £^. Put f = ^jyjdjUj G Ti. Then 
Bf = Y.] djUj = u, and so u G Ran(B) — D(A), as desired. 

Selfadjointness, and discreteness of the spectrum 
Theorem 19.1. A is self adjoint, with domain 

D(A) = Ran(B) = { ^yrVjUj : {q} G f}. 

) 

Furthermore, spec(A) = spec^|j,j,(A) = {yj : j > 1}. 

Proof. We have shown above that A is symmetric and D(A*) = D(A), which 
together imply that A is selfadjoint. 
We will show that if 

AG C\{yi,y2,y3,...} 

then A — A is invertible, so that A belongs to the resolvent set. Thus the spec- 
trum consists of precisely the eigenvalues y,. Note each eigenvalue has finite 
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multiplicity by Theorem 4.1, and is isolated from the rest of the spectrum; 
hence A has purely discrete spectrum. 

The inverse of A— A can be defined explicitly, as follows. Define a bounded 
operator C :7i ^> 7i on f = ^. Cjitj G H by 



]"•)> 



where we note that (yj — A)^ is bounded for all j, and in fact approaches 
as j —} oo, because ly, — A| is never zero and tends to oo as j ^ cxd. This 
new operator has range Ran(C) — Ran(B), because (y, — A)^^ is comparab 



to y= (referring here to the characterization of Ran(B) in Theorem 19.1) 



ThusRan(C) = D(A). 

Clearly (A — A)Cf = f by definition of A, and so Ran(A — A) = %. 
Similarly one finds that C(A — A)u = u for all u G D(A). Thus C is the 
inverse operator of A — A. Because C is bounded on all of "H we conclude 



A — A is invertible, according to the definition in Chapter 17, D 



Example: Laplacian on a bounded domain 

To animate the preceding theory, let us consider the Laplacian on a bounded 
domain O C M'^, with Dirichlet boundary conditions. We work with the 
Hilbert spaces 

and the sesquilinear form 

r 



Vu ■ Vv dx + 
a 



uvdx = (u,v)hi 



Q 



which in Chapter p^ gave eigenfunctions satisfying (—A + 1 )u = (A + 1 )u 
weakly. In this setting, u = Bf means that (—A + l)u = f weakly. Note 
B : \}[D.) -^ Hj(a), and recall that A = B^^ 

Proposition 19.2. The domain of the operator h contains H^(0)nHQ(0), 
and 

A = -A + 1 

onH2(a)nHj(a). 

Furthermore, if 90 is smooth then D(A) = H^[Q) fl Hj(0), in which 
case A = —A + 1 on all of its domain. 
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Proof. For all u G H^(0) n Hj(a),v G Hj(a), we have 

(u,v)hi = (— Au + u,v)i_2 by parts 

= (B(— Au + u),v)hi by definition of B. 

Since both u and B(— Au + u) belong to Hj(0), and v G Hj(Q) is arbitrary, 
we conclude from above that u = B(— Au + u). Therefore u G Ran(B) = 
D(A), and so H^{C1] f] Hj(a) C D(A). 

Further, we find Au — — Au + u because B = A^^ , and so 

A = -A + i on H^(a)nHi(a). 

Finally we note that if 911 is C^-smooth then by elliptic regularity the 
weak solution u of (—A + l)u = f belongs to H^(0), so that Ran(B) C 
H2(a)nHj(a). Thus 

D(A) = H^(a)nHj(a) 

when 90 is smooth enough. In that case A = — A+ 1 on all of its domain. D 
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